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PREFACE 


Objective : 

1. To determine whether ERTS data can be used to characterize 
parameters affecting watershed runoff. 

2. To compare the performance of hydrologic models when 
routine manually determined parameters are used versus when ERTS 
derived parameters are used. 

Scope : 

i 

The study was based on the hydrologic data pertaining to 20 
highly instrumented watersheds located in central Oklahoma. Data 
v from one group of 10 watersheds was related to linear combinations 
of mean digital MSS data to develop a prediction scheme for water- 
shed runoff coefficients. Two storm runoff equations were fitted 
to the watershed data to arrive at measured coefficients that 
represented watershed surface conditions. The SCS storm runoff 
equation was used in this study to illustrate that with dry surface 
conditions the coefficient commonly called the runoff curve number 
can be related to ERTS-MSS digital data. Predictions based on the 
relationship found in two ERTS scenes were tested on the remaining 
10 watersheds. 

Conclusions : 

Predictions were significantly improved over the runoff curve 
numbers calculated by conventional means, and major improvement in 
estimates of flow into flood control works is possible. The use of 
the technique produces objective estimates of watershed runoff that 
can be repeated. This study was not extensive enough to determine 
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if dense vegetation on watershed surfaces will limit the application 
of the technique. 

Summary of Recommendations : 

A follow-on study to test the prediction technique on heavily 
vegetated watersheds and determine the extent of the regions where 
the technique can be used. 
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1. 0 Introduction 

Engineering works for the control of flood waters have 
been built throughout several centuries, but never at the 
pace that will be required to protect life and food production 
capacity for our fast growing world population* Safe and 
economical design of flood control structures such as dams, 
diversion dikes, and conveyance channels require reliable 
estimates of flood flows that can be expected from the water- 
shed areas above the structure. To gain some understanding of 
water supplies and watershed runoff, many hydrologic measure- 
ments; river stages, flow volumes, rainfall and snowfall, and 
topographic measurements have been accumulated in recent times. 
The hydrologic measurements available today pertain primarily 
to large river basins and, in some instances, drainage areas 
of smaller watersheds in regions that have historically high 
density populations. 

Long periods of record, 30 years or more, are desirable 
to establish the rainfall-runoff relations that can be expected 
in both high and low rainfall seasons. Instrumentation and 
manpower to gather and compile adequate hydrologic records are 
expensive and time consuming, therefore , watersheds selected 
for monitoring are usually those where the need for data is 
most critical. The more developed countries of Europe and 
North America have never had resources available to accumulate 
watershed runoff data to meet future needs, while under- 
developed nations have even less data to use in the design of 
flood control structures. 
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The United States Geologic Survey (USGS) has been the 
primary agency responsible for gathering watershed runoff 
information in this country. Their efforts prior to the last 
decade were directed toward data collection for large basins 
where the Department of Interior or the U. S. Army Corps of 
Engineers would be responsible for designing flood control 
works. Runoff from tributary watersheds has been recorded 
by the USGS on only a small percentage of the tributary water- 
sheds as funds would permit. 

The United States Department of Agriculture (USDA) on 
the other hand, has monitored rainfall and runoff from many 
small drainage areas. The small drainage areas usually con- 
sisted of a so-called unit source area having uniform soil and 
cover. Their interest was in the drainage area rainfall-runoff 
relations needed to design small detention and control structures 
for farmers either to control erosion or provide water for live- 
stock. The United States therefore has rather extensive records 
of runoff from most large river basins and numerous small unit 
source areas. Runoff records for agricultural watersheds with 
1.0 to 500 square kilometers in drainage area are generally not 
available . 

Runoff and flood flows from small agricultural watersheds 
have, in recent years, become an important concern of the agri- 
cultural and small town property owners. The control of floods 
by construction of numerous small detention structures on 
branches of a tributary watershed has become an accepted prac- 
tice. At the same time, runoff from small watersheds has become 
important as a water supply for municipalities. 
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Reliable projections of the quantity and rate, of runoff 
from the surface of the land into rivers and streams are 
difficult to obtain for ungaged watersheds. However, this 
information is needed in the design of any structure located 
in the vicinity of a water course; for example, the storage 
capacity of municipal, water supplies and flood detention 
structures. When projections of runoff are questionable, the 
storage capacity of such structures is quite often overde- 
signed. Overdesign not only increases construction costs, 
but may also lead to significant reduction in the flushing 
action needed to maintain good water quality in structures 
where inflow is initially saline and evaporation rates are 
high. 

Examples of overdesign are frequently observed through- 
out the Great Plains. In the central reach of the Washita 
River basin, a water supply reservoir for the city of Chicka- 
sha was completed in 1958. The history of runoff into this 
structure for the past 15 years shows that the anticipated 
runoff exceeded the observed by several times. As a result 
of the overdesign and other related problems, the salinity of 
the stored water has increased sufficiently to preclude use 
as a municipal water supply and even sometimes as an irrigation 
source. 

In an adjacent watershed, Sugar Creek, a study of the 
response of a large number of flood detention structures to 
a large storm on September 19-20, 1965 shows that the inflow 
was only about one-half that expected. 
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The quantities and rates of runoff that are used to design 
structures such as these are estimated from various types of 
watershed runoff models* These models are mathematical equa- 
tions, generally of an empirical form, based on the drainage 
area, topography, soil and cover of the subject watershed. 
Usually, parameters of these models represent drainage area 
and surface conditions prevailing at the, time of a rainfall 
event. The integrated influence of these characteristics com- 
bined with measurable climatological parameters (rainfall, 
antecedent moisture, rainfall intensity! etc.) can produce 
reasonable estimates of storm runoff. Efforts have been made 
to quantify the integrated influence of soils, cover and sur- 
face roughness. However, at the present time, no objective 
means exists for estimating this influence. Present methods 
for estimating this influence are tedious, expensive and 
subject to judgment of the hydrologist. 

The problem of subjectivity associated with present 
methods that are used to derive coefficients of runoff equa- 
tions may be circumvented by applying digitized data obtained 
from ERTS. Even though hydrologic analysis based on photo- 
graphic data obtained from ERTS suffers from the same 
subjectivity as do the existing methods, data obtained from 
the multispectral scanner, MSS, on board the satellite can 
provide digital data that can be incorporated into mathe- 
matical equations. This type data is less subjectively 
biased than are the other types of data. 
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The HRTS multispectral data is also superior to aerial 
digital data collected prior to ERTS because the relatively 
constant sun angle within an ERTS frame,, caused by* the 
satellite's sun synchronous orbit, eliminates many of the 
problems in comparing the spectral reflectance from watershed 
surfaces. The ERTS data has additional advantages in that 
there is sufficient spatial resolution to allow the identifi- 
cation of watersheds as small as 10 hectares, and yet the 
number of data points depicting a large watershed are reduced 
to a manageable number. 

This study was conducted to investigate the possibility 
of using ERTS-MSS digital data to define the coefficients for 
watershed runoff models used on small urigaged watersheds and 
thus provide economical and timely data for planning and 
design of flood detention structures. 

Estimates of runoff from small watersheds are usually 
made either by use of a simple empirical runoff model or by 
use of the modern complex watershed models. The empirical 
watershed models are widely used by practicing hydrologists. 
Numerous brief models have been developed in the last century, 
however, one in particular, developed by the USDA Soil Conser- 
vation Service (SCS) is most commonly used in the United 
States (Mockus, et al., 1971). 

The SCS equation: 


Q - 



CD 
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where S 


1000 


10 


Q * storm runoff (cm/2.54) 

P ■ weighted storm rainfall (cm/2.54) 

S * storage in the watershed surface (cm/2.54) 

CN a function of soil, cover, antecedent moisture 
(dimensionless) 

All of the many empirical runoff equations represent the 
influence of near-surface storage by use of one or more 
coefficients . In this equation, the coeff icient * CN, commonly 
called the curve number is a function of the surface character- 
istics of the watershed at the time a storm event begins. No 
attempt is made to describe the dynamic changes in storage 
through the storm period due to changes in infiltration rate 
in the soil, depletion of available storage in the vegetation, 
or depletion of storage in the soil. 

The SCS runoff equation was developed from data collected 
on small plots located in several regions with a variety of 
soils and cover. Tables were developed from these data to 
define the influence of soil type, vegetative cover, and ante- 
cedent moisture conditions. To apply the equation, each unit 
of area with a single soil type and single cover must be assigned 
its own curve number. The areas and curve numbers are tabulated 
for an entire watershed drainage area and the curve numbers are 
weighted by area to determine an average curve number for the 
watershed. Thus, the computation of the watershed curve number 
is not only dependent on the judgment and experience of the 
hydrologist , but also is time consuming and sub j ect to computa- 
tional errors. 
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The goal of this study has been directed toward develop 
ment of a technique where the coefficient CN for the SCS 
runoff equation can be determined objectively from ERTS data 
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2.0 Approach: 

The problem was approached by first assuming that differ- 
ences in the soil-cover complex over a watershed area would be 
detectable by differences in reflectance of visible and/or 
near infrared light. The reflectance i;(i each band of light 
available from the digital data of the ERTS multispectral 
scanner could then be averaged over a watershed drainage area 
to provide a single value for comparison to measured watershed 
runoff coefficients. 

Twenty watersheds with extensive measurements were 
selected to represent the widest range of rainfall-runoff 
response experienced in central Oklahoma. The watersheds 
were then to be divided into two comparable groups of 10, each 
group having as near as possible, the same range of sizes and 
the same range of runoff coefficients. Group I watersheds 
would then be used to distinguish and develop a relationship 
between the coefficients of a runoff equation and MSS data. 
Group II watersheds were set aside for verification. 

To accomplish this, storm rainfall, storm runoff, rainfall 
intensity, and antecedent rainfall would be calculated and com- 
piled for the available period of record on each of the 20 
watersheds. These data would provide the basis for evaluating 
the actual rainfall-runoff response of each watershed and ; 
determine the measured curve number (CN) for each watershed. 
These data would provide the basis for grouping the watersheds 
and also be sufficient to evaluate other simple empirical 
runoff equations by fitting the data to the equations and 
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optimizing the coefficients. Existing computer programs for 
optimization were available for processing these data. 

Processing of the MSS digital data was planned in several 
steps. First, to elininate the costs of handling large volumes 
of digital data, 70 mm black and white photographs of MSS- 5 
data were requested to screen the available supply and limit 
the study to cloud-free scenes representing all seasons of a 
year. Secondly, computer programs would be developed to identi- 
fy, extract, and isolate the digital data that represented each 
watershed drainage area. A technique would be developed where 
watershed boundaries could be mapped on an overlay for a dis- 
play system allowing selection of coordinates for a series of 
points that would define the boundaries of the location of data 
on the digital tape. Computer programs could then be written 
to excerpt the pertinent data, store it in separate files and 
compute the mean and standard deviation of each band over the 
surface of the individual watershed. 

It was then proposed that the linear combinations of mean 
spectral response from the four bands would be examined by 
multivariate analysis techniques and simple curve-fitting tech- 
niques to find the best relation between the MSS data and 
watershed runoff coefficients. Only data from Group I water- 
sheds would be used to develop the relations. If an acceptable 
relation existed, it was then to be used as a prediction scheme 
on the Group II watersheds. Both predicted coefficients and 
coefficients developed by the conventional SCS procedure would 
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be compared to measured coefficients to determine if the 
remote sensing technique could determine coefficients as well 
or better than the conventional method* 
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3.0 Ground Truth and Data Processing 

3.1 Basic Data 

Ground truth for this study consisted of recorded hydro- 
logic data collected by the Agricultural Research Service 
CARS) from the 1961 through 1972 time period. The 20 water- 
sheds used are located in Grady and Caddo Counties in central 
Oklahoma. The watersheds generally represent small tributary 
watersheds of the Southern Great Plains area. Mean annual 
rainfall in the area is 78 cm. 

Two hundred and fifty-six storm events were selected 
from the records of the 20 watersheds. Storm events were 
selected on a basis of weighted mean storm rainfall greater 
than 3 cm. and measured runoff greater than .03 cm. The 
number of acceptable events ranged from 9 to 21 events per 
watershed. 

3.2 Data Compilation 

Data compiled for each storm event used included weighted 
mean rainfall, runoff, antecedent rainfall index [30-day, 
decayed), antecedent rainfall index [5-day sum), and maximum 
hourly intensity. Drainage area above farm ponds varied from 
0 to 40 percent of the total drainage area within each water- 
shed. Farm ponds would modify runoff to a different extent 
on each watershed, therefore runoff was adjusted to an esti- 
mate of the contributing area using records of farm pond 
storage. The 30-day antecedent rainfall index (Linsley, et al . , 
1949) was computed by depleting the residual rainfall index 
daily by a seasonally varying constant. The constants used 
were derived from an inverse mean daily temperature. 
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Antecedent rainfall index calculated in this manner is usually 
more realistic than simple summation of prior rainfall over 
some period of time. An antecedent precipitation index was 
calculated for the entire study area and values were selected 
for each day an ERTS scene was taken. Antecedent precipita- 
tion index values for the day of each scene used in this study 
are listed in table 1. Summation of rainfall for the 5-day 
period prior to each storm was compiled for these storms since 
the SCS procedure uses this index to account for prior rainfall. 


Scene Number 

ft — — 

Date 

30-Day API (cm/2.54) 

1058 

09-19-72 

.028 

1094 

10-25-72 

2.27 

1184 

01-23-73 

1.94 

1256 

04-05-73 : 

1.19 

1274 

04-23-73 

.928 

1400 

08-27-73 

.0180 

1508 

12-13-73 

1.063 


3.3 Curve Numbers Calculated 

The rainfall and runoff values were used in the SCS runoff 
equation (Equation 1) to calculate actual curve numbers for 
each storm event. It is apparent from a study of these events 
that conversion from one class to another in the SCS routine 
is not appropriate to storms in this study area unless a large 
number of storms have occurred in each antecedent condition 
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class. A large majority of the events were in the Class I 
category of antecedent precipitation index used by SCS. 
Therefore, only Class I storms were used to derive mean curve 
numbers for watersheds in this study. 

3.4 Curve Numbers by Conventional SCS Method 

The Soil Conservation Service had previously determined 
conventional runoff curve numbers for 1 t of the watersheds used 
in this study. These were furnished to ARS along with soils 
maps and photo mosaics of each watershed for computation of 
curve numbers for the remaining 8 watersheds by the conventional 
SCS method. Land use was interpreted from color and color IR 
photographs taken on the ERTS aircraft support flights. Singu- 
lar soils - land use classes were identified, assigned a curve 
number and the area of each class measured. Weighted mean 
curve numbers were then calculated from these data. 

A listing of the curve numbers calculated from measure- 
ments and those calculated by conventional SCS techniques can 
be found in table 2. 

3.5 A Second Runoff Equation 

Attempts were made to fit another runoff equation to the 
data using precipitation, 30-day antecedent precipitation, and 
maximum hourly intensity as variables. Very poor results were 
obtained after trying several linear combinations of the vari- 
ables. Ultimately the intensity was deleted and runoff was 
fitted to rainfall alone, then deviations in predicted runoff 
from measured runoff were fitted to the 30-day decayed 
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Table 2. 

Summary of 

Watershed Data* 




Watershed 

Number 

Drainage 
Area (km^) 

Percent DA 
above Ponds 

Measured 
CN CEq.l) 

Conventional 
CN (Eq.l) 

Constant 
C (Eq.2) 



Group 

1 



206 

.110 

0.0 

53.6 

61 

.034 

207 

.0777 

0.0 

75.8 

86 

.122 

111 

67.3 

26. 4 

60.9 

71 

.038 

141 

190. 

20.0 

58.0 

74 

.023 

512 

91.2 

31.4 

67.2 

74 

.050 

513 

49.7 

34.4 

65.7 

74 

.054 

5141 

16.4 

28.2 

61.5 

74 

.041 

5146 

3.08 

31.1 

63.8 

73 

.068 

522 

539. 

19.5 

57.1 

73 

.031 

612 

2.28 

20.7 

66.7 

74 

.057 



Group 




205 

.0959 

0.0 

54.4 

61 

. 039 

208 

.0749 

0.0 

77.4 

83 

.147 

121 

534. 

21.2 

58.6 

78 

.023 

311 

65.5 

40.7 

69.6 

77 

.078 

511 

154. 

34.4 

69.4 

75 

.082 

5142 

1.39 

45.4 

59.4 

76 

.027 

5143 

1.97 

33.7 

56.3 

68 

.021 

5144 

5.90 

38.4 

62.8 

76 

.066 

611 

19.6 

31.3 

70.2 

77 

.065 

621 

86.2 

20.6 

67.4 

77 

.057 


*A11 values other than drainage area are dimensionless. 


14 


antecedent precipitation thus leading to the exponential values 
to relate runoff to precipitation and antecedent precipitation. 
These exponents were derived using all 256 storm runoff events. 
The resulting equation was in the following form: 

q ;s C p2.15 api . 278 (2) 

in which Q « watershed storm runoff (cm/2.54) 

C ■ a dimensionless coefficient representing 
differences in watershed conditions 
P = weighted mean storm rainfall (cm/2.54) 

API - 30-day decayed antecedent rainfall index 

derived using inverse temperature curves to 
adjust for seasonal variations (cm/2.54) 

The exponents were then fixed in Equation 2 and a mean coeffici 
ent fitted for each watershed (Table 2). The coefficients and 
exponents accepted for this simple equation predict runoff that 
has a multiple correlation with the measured runoff of .7220, 
whereas the SCS equation using curve numbers accepted for this 
study produce a multiple correlation of .7112 when compared to 
the measured runoff. This indicates the two equations used 
are of comparable quality for predicting storm runoff in this 
region. Use of only one or two storm parameters cannot be 
expected to produce better results than this. 

3.6 Map Requirements 

Maps were obtained for each of the watershed areas. The 
coordinates of a series of points defining the watershed 
boundaries and major stream channels were selected on a chart 
reader to produce a card deck for each watershed. These data 
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decks were stored on disk files and used as control for a 
plotter program to produce overlay maps. Overlays were then 
plotted to match the scale of ERTS data displayed on a tele- 
vision screen or to match a conventional grey scale computer 
printout. Location of major water bodies in or near the 
watersheds were also mapped to aid in positioning the boundary 
overlay. 
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4.0 ERTS Data Processing 

4.1 Data Screening 

Microfilm of the ERTS scenes and the data search system 
available at Goddard Space Flight Center were used for pre- 
liminary screening to select relatively cloud-free scenes 
over the study area. Nine-inch photographs of the data were 
then used for a second look before bulk digital tapes of the 
MSS data were ordered, 

4.2 Digital Data Selected 

Mul tispectral digital data for the watershed areas were 
obtained from the ERTS scanner for seven scenes that covered 
the study area. Each scene represented major changes in soil 
moisture and vegetative conditions. The first scene, 1058, 
dated September 19, 1972, shows a dry dormant condition with 
almost no ground cover. The second scene, 1094, dated October 
25, 1972, provides data with essentially the same ground cover 
but extremely wet conditions. The third scene, dated January 
23, 1973, showed minor growth in winter wheat fields and wet 
conditions. The scenes, 1256, dated April 5, 1973 and 1274, 
April 23, 1973, showed substantial ground cover for crops, 
pastures, and timber. One of the growing season scenes was 
extremely wet and one moderately dry. The last two scenes, 
dated August 27, 1973 and December 13, 1973 show extremely 
heavy vegetative growth on cropland for the fall and winter 
seasons. Table 1 lists the calculated 30-day antecedent 
precipitation index associated with each scene used. Reflec- 
tance in the near infrared, band 7, was relatively high in 
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the grassland areas on the fall scene, indicating more growing 
•vegetation than on the scene from the previous fall . 

4.3 Data Processing Approach 

The MSS data from ERTS were obtained in the form of 
sequential and adjoining tapes that were laced together with 
the aid of a computer program (MERGE, Appendix) forming a 
single file for each watershed area. Due to the nature of 
the problem, it is very important to keep the relative position 

of the data points correct and to be able to accurately define 

its location within a watershed. Without this, coordination 
of ground truth and spectral response would be impossible. 

Since the methods used to display the MSS digital data present 
an image enlarged in the cross -track axis, the maps developed 
as part of the ground truth data were also enlarged in the 
cross-track axis before they were used to locate watershed 
boundaries. 

Data sets for the larger watersheds were obtained by dis- 
playing the entire ERTS frame for the area of interest on a 

television display. The portion representing the watershed 
was isolated using the distorted maps and a computer program 
(OKLAH, Appendix). The computer program represents the water- 
shed boundary by a series of adjacent parallelograms. About 
20 parallelograms seem to be enough to adequately define the 
watershed boundaries (Fig. 1). This program stores the 
watershed data on a secondary tape in a rectangular file 
placing zero values in all data points outside the watershed 
boundary. The ERTS data format is retained on the secondary 
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tape so the data pertaining to the watershed areas can be 
displayed on the screen for visual verification that the data 
selected covers the entire watershed drainage area. 

Data sets for the small watersheds were more easily 
obtained by displaying the ERTS data on computer printout and 
overlaying it with distorted maps. In both procedures, ponds 
and channels aid significantly in locating specific points 
that assure proper selection of the watershed boundary points 
on the MSS data tapes. In some data sets, stream channels are 
difficult to locate in grassland areas; a study of the ERTS 
data showed that MSS band 5 resolved the stream channels and 
ponds sufficiently well to position the overlay. However, 
during the growing season, (MSS band 4 + MSS band 5)/MSS band 6 
enhanced the resolution of the stream channels. Although 
cumbersome, the system described for small watersheds will work 
to locate larger watersheds if a display system is not available. 
Ratioing bands 5 and 7 also helped to enhance some scenes where 
single bands offer little contrast. 

A simple computer program (Mean 4, Appendix) was used to 
calculate the mean and standard deviation of the digital values 
for each band from the digital values stored on the secondary 
tapes. Computations were made for each watershed on all scenes 
used in the study. These values were considered as the basic 
set of multispectral scanner data that would be compared to 
the hydrologic data. A summary of these data from the seven 
scenes can be found in the Appendix. 

The MSS data from scene 1058 was used as a base to deter- 
mine if all data points were necessary in the computation of 
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mean values for each watershed. From each of the watersheds 
represented by more than 5,000 data points, 256 independent 
samples were selected. By calculating the means, sample size, 
and standard deviation of each sample, then combining adjacent 
pairs of samples and repeating the process, the change in the 
mean and standard deviation with change in sampling frequency 
was observed. Less than 1 percent change in the mean occurs 
between any sample greater than approximately 2,500 data points 
and the sample that includes all data points in the watershed. 

In an operational system it seems reduced sampling could there- 
fore be used to cut computer costs when studying large watersheds 
without impairing the quality of the data. 

4.4 Aircraft MSS Data 

Aircraft tapes from the 24-channel MSS flown in support of 
this project were converted to an ERTS-type format for display 
on the Dicomed. The data quality seems to be erratic and 
unusable for a complete watershed area. 

The aircraft MSS data was, however, useful for a data base 
in an incidental water quality study. Digital values were 
punched on cards for a few selected areas that contained ponds 
where water quality samples were collected near the time of 
the flight. A technique for processing of the data over water 
bodies was developed where distribution of the digital values 
in each band was defined and only the values falling within 
one standard deviation from the mean were used in the analysis. 
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5.0 Analysis 


5.1 Multiple Discriminant Analysis 

Since the objective of this study is to see if MSS data 
can be used to calculate parameters of a runoff equation, 
several means of relating the two data sets were used. Dis- 
criminant analysis, which can be used to study group similarity 
or difference and relate this to group descriptors, was used 
to examine the MSS data of the watersheds having extreme 
differences in observed runoff coefficients. Using a modified 
multiple discriminant analysis program [Cooley and Lohnes, 1962), 
good discrimination was observed between these watersheds when 
MSS bands 4, 5, and 7 were used in the linear discriminant 
function. The good discrimination was found in the dry dormant 
scenes, however the same band also produced the best discrim- 
nation in other scenes. 

Multiple discriminant analyses considering each of the 10 
developmental watersheds as an independent group showed very 
significant group discrimination. However, the discrimination 
did not appear to be related to the runoff coefficients. There 
was no multiple discriminant analyses program available that 
would evaluate two groups of 10 watersheds and maintain the 
ranking of the dependent variable (in this case, the runoff) 
within each group. 

5 . 2 Curve Fitting 

Alternatively, plots of the mean value for each band vs. 
the observed runoff coefficient were made. Since it appeared 
there might be a relationship, all possible combinations of 
the means were plotted. Two promising relationships were 
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evident. The mean of MSS band 5 (u s ) minus the mean of MSS 
band 4 (u 4 ) is reasonably well related to the SCS runoff curve 
number for the watersheds in Group I. This relationship is 
most evident in the dry dormant scene 1058. The combination 
u 5 + y 6 “ ^ u 4 + Produces a more consistent relationship 

when all scenes are considered. Both linear combinations of 
the bands were used to define prediction curves for scene 1058 
CFig. 2). 

Scene 1058 occurred at a time when the 30-day antecedent 
precipitation index was extremely low. Data from scene 1400 
acquired nearly a year later also represented extremely dry 
conditions. Live vegetation was more prevalent in scene 1400 
than in the prior year, however, the 30 -day antecedent rainfall 
index was .028 and .018 for scenes 1058 and 1400, respectively. 
Some cloud cover was found on scene 1400 and watersheds 611 and 
612 could not be identified. Therefore, data from only nine 
watersheds in each group are available for analysis of this 
scene. 

Figure 3 illustrates that a similar relation to the one 
found in scene 1058 can be described by using eight of the nine 
data points. The one data point that plots as an outlier 
belongs to watershed 111. The shift in this point is likely 
due to the influence of cloud cover, therefore it was not con- 
sidered when locating the curves. A summary of the data used 
to plot the prediction curves is presented in table 3. 
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MEASURED CURVE NUMBER 


Figure 2. The relations of MSS data from Scene 1058 
to measured watershed runoff curve numbers. 
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Figure 


54 58 62 66 70 74 

3. The relations of MSS data from Scene 1400 to 
measured watershed runoff curve numbers. 
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Table 3. Data 

Base 

for Prediction 

Curves - 

Group 

I 






Watershed No. 


205 

207 

111 

141 

512 

513 

5141 

5146 

522 

612 

Scene 1058 

HS-M4 


.18 

5.17 

.49 

.43 

1.66 

1.09 

.81 

.75 

-.11 

1.97 

w S +, 1 6' , V 2u 7 


.60 

7.18 

.74 

.18 

2.36 

1.56 

1.34 

1.28 

.07 

2.96 

Coefficient 1 
(Measured CN) 


54.4 

75.8 

60.9 

58.0 

67.2 

65.7 

61.5 

63.8 

57.1 

66.7 

Coefficient 2 
(Conventional 

CN) 

61 

86 

71 

74 

74 

74 

74 

73 

73 

74 

Scene 1400 

u 5 -w 4 


-5.80 

.00 

-2.36 

-4.68 

-3.95 

-3.92 

-4.72 

-4.91 

-5.38 


w 5 +u 6 -W4- 2 w 7 


-7.62 

1.32 

-1.98 

-5.47 

-3.85 

-3.75 

-4.71 

-5.71 

-6.16 




5.3 Prediction of Curve Numbers for Test Watersheds 

Using data from scene 1058 , both relationships illustrated 
in figure 2 were verified on the Group II watersheds (Figs. 4 
and 5). In figure 4, using two bands of MSS data, predictions 
deviated an average of 4.13 units (absolute) from the measured 
values. In figure 5 when the predictions were based on using 
four bands of data, they deviated an average of 3.17 units 
(absolute) from the measured values. 

Curve numbers were then predicted for the Group II water- 
sheds by using the relationships developed from data for 
scene 1400. The predicted values and the conventional SCS 
curve numbers were plotted versus the measured curve numbers 
(Figs. 6 and 7). The average deviation of the predicted values 
from the measured curve numbers was 4.59 units (absolute) when 
using two bands of MSS data and 3.70 units (absolute) when using 
four bands of data. The average deviation of the predicted 
values can be compared to an average deviation of 10.72 between 
the conventional and measured values. 

The predicted curve numbers for the Group II test watersheds 
are summarized in table 4. These data were used to plot figures 
4 , 5 , 6 and 7 . 

Similar plotting techniques were used to examine the 
relation of the coefficient for equation 2. The results were 
comparable, however the relationships are not well defined and 
none seem as promising as the relationships found using 
equation 1. 
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CONVENTIONAL AND PREDICTED CURVE NUMBERS 






Figure 5. Comparison of conventional SCS curve numbers 
to curve numbers predicted with 4 bands of 
ERTS-MSS data (Scene 1058) 
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CONVENTIONAL AND PREDICTED CURVE NUMBER 



Figure 6. Comparison of conventional SCS curve numbers 
to curve numbers predicted with 2 bands of 
ERTS-MSS data (Scene 1400). 
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CONVENTIONAL AND PREDICTED CURVE NUMBER 


82 h • Conventional curve number 

| O Predicted curve number 
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MEASURED CURVE NUMBER 

Figure 7. Comparison of conventional SCS curve numbers 
to curve numbers predicted with 4 bands of 
ERTS-MSS data (Scene 1400) . 
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Table 4. Measured and Predicted Runoff Coefficients - Group II 

Watershed No. 206 208 121 311 511 5142 5143 5144 611 


Scene 1058 
Predicted CN 


05 -V 4 ) 

49.6 

69.0 

62.2 

72.9 

70.6 

62.9 

49.0 

60.6 

64.6 

69.6 

Predicted CN 

(u 5 +U 6 -U 4 - 2 y 7 ) 

51.6 

68.8 

58.1 

72.8 

69.8 

63.9 

49.1 

62.9 

65.6 

68,5 

Scene 1400 
Predicted CN 
(ys‘^4) 

56.7 

65.8 

67.9 

71.6 

72.8 

53.3 

58.4 

62.9 


68.0 

Predicted CN 

05 +^- 114 - 2 * 17 ) 

55.0 

67.2 

65.1 

75.1 

71.9 

59.8 

61.6 

62.7 

- - 

68.8 

Measured CN 

53.6 

77.4 

58.6 

69.6 

69.4 

59.4 

56.3 

62.8 

70.2 

67.4 

SCS CN 

61. 

83. 

78. 

77. 

75. 

76. 

68. 

76. 

77. 

77. 



5.4 Secondary Testing of the Prediction Scheme 

In the introduction it was mentioned that runoff into a 
large number of flood detention structures on Sugar Creek was 
normally about one-half what was expected when the structures 
were designed. To further check this prediction scheme, the 
digital data for the subwatersheds on Sugar Creek watershed 
No. 121 were examined. Data from scene 1058 was used. A grey 
scale map of the Sugar Creek area was printed and an overlay 
was used (Fig. 8) to select data points within each small water- 
shed. The mean difference between values of band 5 and band 4 
was calculated for each subwatershed. 

Long-term hydrologic data on the subwatersheds was not 
available. A severe storm had occurred over the area in 1965 
where rainfall records from a 3-mile rain gage grid were 
available and high water marks in the flood detention structures 
had been collected. The change in storage volume and the 
weighted rainfall were used calculate the measured runoff 
curve numbers. Curve numbers based on this type of data are 
not considered as reliable as curve numbers calculated from 
long-term records. Using the mean difference calculated for 
each subwatershed from the MSS data as entry points to the 
prediction curve in figure 2A, a predicted runoff curve number 
for each subwatershed was obtained. The conventional SCS curve 
numbers and predicted values were plotted versus the measured 
curve number (Fig. 9) to illustrate the improvement possible 
by using the ERTS, MSS data. 
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Figure 8. Grey scale printout with map overlay of the 
Sugar Creek subwatersheds used to test the 
prediction of runoff curve numbers with 
ERTS-MSS data. 
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CONVENTIONAL AND PREDICTED CURVE NUMBER 



Figure 9. Comparison of SCS design watershed runoff 
curve numbers to curve numbers predicted 
using 2 bands of ERTS-MSS data. 
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The average absolute deviation of predicted values from 
the measured runoff curve numbers was 10.18. Average devi- 
ations of the conventional curve numbers for the subwatersheds 
was 24.08. Figure 9 shows that only 3 of the 22 subwatersheds 
have predicted values slightly under their measured value. 
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6.0 Comments and Conclusions 

At the time this experiment was proposed, several factors 
important to the outcome were not fully recognized* First, it 
seemed that the multiple discriminant analysis computer programs 
could be used to find the ultimate equation that would define the 
relationship of ERTS-MSS data to watershed runoff coefficients. 

These programs will identify the variables that obtain best dis- 
crimination between members of a group, however, the dependent 
variable is not confined to a ranking. When more than two sets 
of watershed data are used, a linear combination of the independent 
variables may produce the best discrimination, but result in a 
change in the ranking of the dependent variable. The system, 
therefore, is excellent as a search routine using the high and low 
runoff -producing watersheds to indicate feasible linear combinat ions 
of the MSS data. In its present form the program does not necessarily 
indicate the best linear combination for prediction of the dependent 
variable when several watersheds are considered. 

Secondly, it was assumed that an empirical equation including 
more storm parameters than rainfall amounts would predict storm 
runoff better than the SCS runoff equation. The data available for 
the 20 watersheds was not adequate to develop an equation that pro- 
duced improved results. Equation 2 includes the influence of 
antecedent precipitation which is not incorporated in the SCS 
equation, and still does not produce better estimates of storm 
runoff than the SCS equation. 

The quality and number of usable scenes of the ERTS-MSS data 
was also underestimated prior to the launch. In fact, it now seems 
that the lack of adequate long-term rainfall and runoff records 
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on small watersheds may, in some instances, be the limiting factor 
in making full use of MSS data for watershed runoff estimates. To 
calibrate the ERTS data and apply the technique developed in this 
study, records are necessary on several watersheds with a range of 
runoff producing capability from the lowest to highest runoff curve 
numbers within each ERTS scene. 

The similarity of the curves in figures 3 and 4 indicate that 
the technique used in this study is repeatable in scenes where dry 
surface conditions exist. The slope of the curves in these figures 
also illustrate the fact that the prediction curves based on spec- 
tral differences using all four bands of data are more sensitive 
than curves based on two bands. The digital values in all bands 
for scene 1400 were very large in comparison to the digital values 
in scene 1058. The shift of digital values from one scene to 
another is evident in all seven scenes, but the shift does not 
seem related to surface moisture conditions. 

The good relationship between the linear combinations of digital 
data occurred only in the dry scenes. Table 1 shows that scenes 
from the late spring or summer when vegetative growth was heavy 
were all representative of relatively wet conditions. These data 
are inadequate to isolate the influence of vegetative growth, thus, 
heavy vegetation must at present be considered as a possible limi- 
ting factor in application of this technique. These data do 
indicate that wet surface conditions limit the application of the 
technique during the dormant season. Scenes 1094 and 1184 repre- 
sent essentially the same sparse vegetative cover that existed in 
scene 1058, but both winter scenes represent wet surface soils. 
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The relationships developed with scene 1058 data do not exist in 
the wet scenes. 

When the prediction scheme was tested on, the Group II water- 
sheds, both linear combinations of MSS bands predicted runoff curve 
numbers better than the conventional calculated curve numbers. 

When predictions based on two bands of data (Figs. 4 and 6) are 
compared to predictions from all four bands (Figs. 5 and 7), it is 
apparent that the four-band system underpredicts less frequently. 
Underprediction of curve numbers by more than five units may pro- 
duce unexpected high flows through the emergency spillway of flood 
detention dams that may peopardize the structure. The most desirable 
prediction scheme would produce values that would fall on or 
immediately above the optimum design line in these figures. To 
reduce the risk of underdesign in an operational system, the pre- 
diction curve (Figs. 3 and 4) could be shifted to the right until 
all data points fall on the left side of the curve. 

The testing of the predictions for the 22 subwatersheds of 
Sugar Creek offers more substantial proof of the validity of the 
technique. These subwatersheds are all above existing SCS flood 
detention reservoirs and represent a range of drainage area size 
common to structures built by SCS under present laws. Figure 9 
illustrates that predictions based on ERTS data reduced overesti- 
mation of curve numbers by more than a factor of 2.36. Curve 
numbers underpredicted were 1.2, .45, and 3.0 units below curve 
numbers determined by the one severe storm. These underpredictions 
would not be considered serious. 

Due to the form of the SCS equation, the influence of the 
improved curve numbers is nonlinear. If a watershed has an actual 
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curve number of 50, an overestimation of 20 units will produce a 
calculated runoff of approximately 7.12 cm from a 15.24 cm rain- 
fall, while an overestimation of 10 will produce calculated runoff 
of approximately 4.8 cm. The actual runoff would be approximately 
2.9 cm. The storage volume would be overestimated 146 percent for 
20 units change in curve number and only 68 percent for a change 
of 10 units. The improvement in prediction of runoff is greater 
than the linear improvement in the estimation of the curve numbers. 

In summary, this study has shown that when dry surface con- 
ditions exist, linear combinations of MSS digital data can repeatedly 
be related to the watershed runoff coefficient used in the SCS storm 
runoff equation. Predictions based on the relationship between 
ERTS-MSS data and measured watersheds can improve SCS runoff curve 
numbers by more than a factor of two over curve numbers calculated 
by the subjective conventional methods . The improvement in esti- 
mating curve numbers can significantly improve estimates of runoff 
necessary for the design of flood control structures. The technique 
developed in this study may be limited by reflectance from dense 
vegetation and should be tested to define the influence of dense 
vegetation on the results. 
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APPENDIX 



SUMMARY OF ERTS-MSS MEANS • ST ANDARD DEV I AT I ON t RUNOFF COEFFICIENTS 





ERTS 

ORBIT 1 

DUMBER 1058 < 

GROUP I 

ON 9 

19 72 




ws 

MSS-4 

MSS-5 

MSS-6 

MSS-7 

SD-4 

SD-5 

SD-6 

SD-7 

COEF1 

COEF2 

COEF3 

205 

31.00 

31.18 

36.36 

17.97 

1.64 

2.21 

2.21 

0.77 

54.4 

61. 

0.039 

207 

32.09 

3 7.26 

39.83 

18.91 

1.78 

3.45 

2.53 

1.00 

75.8 

86. 

0.122 

111 

31.39 

31.88 

39.73 

19.74 

4.38 

7.89 

6.30 

3.15 

60.9 

71. 

0.038 

141 

28.85 

29.28 

37.81 

19.03 

3.36 

7.34 

1.01 

3.34 

58.0 

74. 

0.023 

512 

29.72 

31.38 

36.08 

17.69 

2.45 

5.38 

4.01 

1*88 

67.2 

74. 

0.050 

513 

29.16 

30.25 

35.57 

17.55 

2.36 

5.14 

3.66 

1.66 

6 5*7 

74. 

0.054 

5141 

29.43 

30.24 

34.71 

17.09 

2.09 

4.12 

3.00 

1.29 

61.5 

74. 

0.041 

5146 

29.41 

30.16 

35.07 

17.27 

2.20 

4.00 

2.67 

1.11 

63.8 

73. 

0.063 

522 

30.39 

30.28 

37.04 

13.43 

4.80 

7.39 

4.99 

2.32 

57.1 

73. 

0.031 

612 

30.74 

32.71 

38.71 

18.86 

2.39 

GROUP 

4.77 
1 I 

4.76 

2.68 

66*7 

74. 

0.057 

WS 

MSS-4 

MSS-5 

MSS-6 

MSS-7 

SD-4 

SD-5 

SD-6 

SD-7 

COEF1 

COEF2 

COEF3 

206 

30.63 

30.19 

36.19 

18.06 

1.41 

2.04 

1.94 

0.62 

53-6 

61. 

0.034 

208 

30.08 

32.38 

35.81 

17.42 

2.08 

3.14 

2.79 

1.10 

77.4 

83. 

0.147 

121 

30.44 

31.08 

41.04 

20.64 

8.37 

11.21 

10.38 

5.46 

58.6 

78. 

0.023 

311 

32.49 

36.32 

38.66 

18.57 

2.90 

5.31 

4.57 

2.23 

69.6 

77. 

0.078 

511 

30.72 

33.61 

38.52 

18.82 

2.82 

6.02 

4.64 

2.16 

6 9.4 

75* 

0.082 

5142 

29.14 

29.92 

34.80 

17.14 

1.77 

4.20 

3.27 

1.22 

59.4 

76* 

0.027 

5143 

28.86 

28.30 

32.88 

16.51 

1.65 

2.82 

2*27 

1.10 

56.3 

68. 

0.021 

5144 

29.57 

30.28 

34.80 

17.12 

2.22 

4.03 

2.89 

1.24 

62-8 

76. 

0.066 

611 

32.32 

33.40 

38.65 

18.87 

3.50 

5.42 

5.02 

2.70 

70.2 

77 . 

0.065 

621 

WS 

2 9.72 

32.28 38.40 18*93 

WATERSHED NUMBER 

2.68 

6.39 

5.03 

2.33 

67.4 

77 . 

0.057 


COEF1 — SCS CURVE NUMBER CALCULATED FROM MEASURED RAINFALL AND RUNOFF 
COEF2 — SCS CURVE NUMBER DERIVED BY CONVENTIONAL SCS TECHNIQUE 
COEF3 — COEFICIENT FOR EQUATION 2 
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SUMMARY OF ERTS-MSS MEANS ♦ S TANQARD DEVIATION tRUNOFF COEFFICIENTS 
ERTS ORBIT NUMBER 1094 ON 10 25 72 
GROUP I 


WS 

MSS-4 

MSS-5 

MSS-6 

MSS-7 

SD-4 

SD-5 

SD-6 

5D-7 

COE F 1 

COEF2 

COEF3 

205 

29*96 

29.13 

30.83 

15.30 

1.52 

1.84 

1.53 

0.88 

54*4 

61* 

0.039 

207 

28.07 

29.07 

29.20 

14.20 

0.88 

2.40 

1.66 

1.C1 

75.8 

86. 

0.122 

111 

27.72 

26.23 

30.13 

15.00 

2.65 

4.92 

4.19 

2.17 

60.9 

71 . 

0.038 

141 

26.75 

25.56 

29.28 

14.76 

2.49 

4.62 

4.15 

2.46 

58.0 

74. 

0.C23 

512 

26.65 

25.45 

27.35 

13.42 

1.62 

3.00 

2*61 

1.43 

67.2 

74. 

0.0 50 

513 

26.45 

25.07 

27.18 

13.37 

1.62 

2.94 

2.52 

1.40 

65*7 

74. 

0.054 

5141 

26*65 

25.15 

26.32 

13.15 

1.48 

2-38 

2.04 

1.12 

61.5 

74. 

0.C41 

5146 

26.33 

24.71 

26.74 

13.16 

1.49 

2.32 

1.91 

1.05 

63.8 

73. 

0 • 0 6 w 

522 

26.93 

24.90 

27.15 

13.32 

2.42 

3.63 

3.82 

2.17 

57-1 

73. 

0.031 

612 

26*88 

25.11 

27.62 

13*76 

1.76 

3*01 

1.94 

1.07 

66*7 

74. 

0 • C 5 7 






GROUP 

I I 






WS 

MSS-4 

MSS-5 

MSS-6 

MSS-7 

SD-4 

SD-5 

SD-6 

SO* -7 

COEF1 

COEF2 

COEF3 

206 

28.45 

26.90 

29.07 

14.66 

1.09 

1.61 

1.93 

0.77 

5 3.6 

61. 

0.034 

208 

26.83 

26.03 

27.67 

13.43 

1.44 

1.79 

1.97 

1.10 

77.4 

83. 

0.147 

121 

26.84 

26.21 

30.08 

15.14 

2.77 

5.58 

4.67 

2.48 

58.6 

78. 

0*023 

311 

27.20 

26.42 

27.10 

13.19 

2.11 

3.28 

3.99 

2.23 

6 9*6 

77. 

0.078 

511 

27.12 

26.57 

28.36 

13.91 

1.80 

3.15 

2.94 

1.57 

69*4 

75. 

0.032 

5142 

26.38 

25.90 

27.34 

1 3 « 41 

1.33 

2.54 

2.30 

1.11 

59.4 

76. 

0.027 

5143 

26.68 

24.88 

26.20 

12.95 

1.39 

1.98 

1.75 

0 . 93 

56.3 

63. 

0.021 

5 144 

26.58 

25*00 

26.31 

13.14 

1.52 

2.39 

1.97 

1. 10 

6 2*6 

76. 

0.066 

611 

26.47 

24.42 

27.32 

13.51 

1.47 

2.68 

2.42 

1.40 

70.2 

77. 

0.065 

621 

26.17 

25.26 

27.91 

13.82 

1.65 

3.45 

2.88 

1.54 

67.4 

77. 

O.Q37 


WS WATERSHED NUMBER 

C0EF1 — SCS CURVE NUMBER CALCULATED FROM MEASURED RAINFALL AND RUNOFF 
C0EF2 — SCS CURVE NUMBER DERIVED BY CONVENTIONAL SCS TECHNIQUE 
COEF3 — COE FI Cl ENT FOR EQUATION 2 



SUMMARY OF ERTS-MSS MEANS » STANDARD DEV I AT I ON « RUNOFF COEFFICIENTS 
ERTS ORBIT NUMBER 1184 ON 1 23 73 
GROUP I 


ws 

MSS-4 

MSS-5 

MSS**6 

MSS«*7 

SD~4 

SD“5 

SD"6 

SD-7 

COEF1 

CQEF2 

CCEF3 

205 

20*50 

20*96 

22.62 

11*88 

1.14 

1.04 

1.92 

1.28 

54.4 

6 1 e 

0.039 

207 

1 9 • 5 2 

21*10 

21*61 

11 • 52 

0.68 

1*58 

1.25 

0.93 

75.8 

8 6 a 

0.122 

111 

21.06 

21.04 

22.77 

11.83 

2.23 

5 « 6 6 

3.99 

2*16 

60*9 

71 * 

0.C38 

141 

19*88 

19*75 

21.55 

11.42 

2*09 

4*05 

4.39 

2.48 

58.0 

74 * 

0.023 

512 

1 9 o 62 

19*14 

20*63 

10.77 

lolQ 

2e 12 

2.93 

1.81 

67.2 

74 * 

0.050 

513 

19*58 

19*08 

20.52 

10*74 

1*11 

2*13 

3*02 

1.89 

65.7 

74 * 

0.054 

5141 

19*64 

19*00 

20.02 

10*41 

1*08 

1*83 

2*13 

1.33 

61*5 

74. 

0*041 

5146 

19 s 57 

18*83 

19.77 

10.35 

1*10 

2*00 

2.21 

1*33 

63.8 

7 3 ® 

0*06-3 

522 

20*66 

20*15 

22.43 

11.84 

2 a 39 

3 o 92 

4 a 2 4 

2 e 30 

57® 1 

73 o 

06031 

612 

1 9 o 69 

1 9 g 4 5 

2 1 o 19 

11.25 

1*17 

2.52 

3 o 2 1 

1 0 86 

6 6 o 7 

74 e 

0*057 






GROUP 

I I 






WS 

MSS = 4 

MSS®*5 

MSS"6 

MSS“7 

SD“4 

SD“5 

SD“6 

SD-7 

COEF1 

COEF2 

COEF3 

206 

20 o 41 

20 o 48 

22.52 

12 a 07 

0.82 

0 s 9 1 

1.57 

1 o 13 

5 3 o 6 

6 1 o 

0 o 0 3 4 

208 

19*31 

1 9 o 00 

19 e 88 

1 0 o 1 9 

0 o 93 

1*36 

1 Q 1 8 

0.98 

7 7 o 4 

8 3 g 

0® 14? 

121 

1 9 o 8 3 

1 9 o 8 6 

21*97 

11.54 

2 o 33 

4o 7 1 

4 o 4 7 

2 o 33 

5 8 « 6 

78s 

0*023 

311 

1 9 o 78 

19 o 50 

22.01 

1 1 o 6 1 

1*18 

2 o 1 6 

2*44 

1 a 46 

69 o 6 

77 s 

0 a 07 8 

511 

1 9 o 8 6 

19a77 

21*85 

1 1 o 45 

1 o 33 

2.32 

2*84 

1.69 

69 a 4 

75 o 

0.082 

5142 

20 a 7 1 

2 0 g 6 5 

22 o 37 

1 1 o 75 

1 o 52 

2 o 2 7 

2 o 48 

1 o 4 5 

5 9 o 4 

76 o 

0.027 

5143 

20 o 1 0 

19 o 8 5 

2 0 a 70 

10.76 

0 o 9 7 

1.35 

1 o 9 7 

1 o 33 

5 6 o 3 

68 o 

0. 021 

5144 

1 9 o 5 1 

1 8 o 7 1 

1 9 o 72 

10o 28 

1*10 

1 o 9 7 

2.18 

1.38 

62 o 8 

7 6 o 

0.066 

611 

1 9 o 73 

1 9 o 2 0 

21 *12 

1 1 o 03 

1.18 

2 ° 1 9 

2 o 9 2 

1.79 

70 * 2 

77 o 

0.065 

621 

1 9 o 5 8 

1 9 o 8 1 

21 o26 

1 1 o 09 

1*18 

2 o 7 5 

3.20 

1 o 83 

6 7 o 4 

77 o 

0.057 


W$ — — WATERSHED NUMBER 

COEF1 = - SCS CURVE NUMBER CALCULATED FROM MEASURED RAINFALL AND RUNOFF 
COEF2 “= SCS CURVE NUMBER DERIVED BY CONVENTIONAL SCS TECHNIQUE 
COEF3 COE F I C I ENT FOR EQUATION 2 



SUMMARY OF ERTS~MSS MEANS ^STANDARD DEV I AT I ON & RUNOFF COEFFICIENTS 
ERTS ORBIT NUMBER 1256 ON A 5 73 


ws 

MSS-4 

MSS“ 5 

MSS“6 

MSS~7 

GROUP 

SD°4 

I 

SD~5 

SD“6 

SD~7 

COEF1 

COEF2 

COEF3 

205 

30.81 

30 e 78 

41*30 

22*81 

1 a 04 

2*10 

2.71 

1 * 44 

5 4 e 4 

61® . 

0.039 

207 

34.38 

43 o 4 3 

49*81 

26*29 

2 c 40 

5*26 

3*17 

1*79 

75*8 

86* 

0.122 

111 

33*18 

34 o 4 1 

45*85 

25.06 

4 o 2 7 

8 * 56 

7 * 5 1 

4*55 

60.9 

71. 

0.038 

141 

3 1 e 3 4 

34*02 

42 o 46 

22 o 90 

3*62 

8*02 

8*57 

4.90 

58*0 

74 o 

0.023 

512 

30*21 

31*03 

38.54 

20*68 

1 e 8 9 

4o 56 

7e05 

4.69 

67 o 2 

74 o 

0.050 

513 

30*18 

30 s 98 

3 7 o 06 

19.72 

1 e 77 

3 o 9 3 

5 * 9 8 

3 s 90 

6 5.7 

74. 

0.054 

5141 

30*18 

30o 72 

35*53 

X 8 e 82 

1 o 42 

2 s 73 

4 o 1 0 

2 s 48 

6 1 e 5 

74. 

0.041 

5146 

30*16 

3 0 o 8 3 

35*61 

18*92 

1*55 

2*70 

4*11 

2*56 

63*8 

73 o 

0*063 

522 

32*65 

32*71 

43 o 59 

23*89 

7*58 

7*65 

7*36 

4*31 

5 7 o 1 

73* 

0*031 

612 

30 o 50 

31 *29 

42*20 

2 2 o 97 

1*73 

4*23 

6*05 

3*95 

6 6*7 

74* 

0*057 

ws 

MSS-4 

MSS“5 

M$S“6 

MSS“7 

GROUP 

SD-4 

I I 

SD m 5 

SD~6 

SD a 7 

COE F 1 

CQEF2 

COEF3 

206 

30*84 

31*06 

39*81 

2 1 o 84 

1*10 

3 o 09 

2*15 

1 o 32 

53*6 

6 1 * 

0*034 

208 

30 o 2 2 

33*22 

39 o 1 9 

20o 63 

1*19 

3*43 

4 o 2 8 

2 o 66 

77*4 

8 3 o 

0 * 14 7 

121 

32*55 

36*40 

45*71 

24*48 

4*53 

9*96 

9*31 

5 s 10 

58*6 

78 * 

0*023 

311 

31*06 

31*22 

4? o 40 

26*36 

3 o 3 5 

7*77 

6*05 

4* 28 

6 5*6 

77 o 

0 o C7S 

511 

30*62 

31*39 

43 o 59 

2 3 o 76 

2*93 

6 o 98 

6*92 

4*67 

69*4 • 

75* 

0.082 

5142 

30 o 2 7 

30*01 

34 a 52 

18*31 

1*33 

2*46 

3*33 

2*62 

59*4 

76* 

0*027 

5143 

30*45 

30*63 

34 o 88 

18o 36 

1*16 

2* 30 

3*45 

2*15 

56*3 

68 o 

0 o 02 1 

5 144 

30*05 

3 0*64 

35*61 

18*87 

1*46 

2 * 8 G 

4 o 2 2 

2*61 

62*8 

76 o 

0.066 

611 

29 o 9 8 

30*43 

3 8 o 38 

20*90 

1*91 

4*66 

5 o 99 

4 o 0 7 

70*2 

77 o 

0 o 065 

621 

30*84 

33 * 94 

4 1 o 1 0 

21*74 

2*15 

5 o 8 7 

6*90 

4*09 

67*4 

77* 

0*057 

tyy! ^ taacseaa 

COEF1 — 

WATERSHED NUMBER 

SCS CURVE NUMBER CALCULATED 

FROM MEASURED 

RAINFALL AND 

RUNOFF 

COEF2 “«= 
C0EF3 

SCS CURVE NUMBER DERIVED BY 
COEFICIENT FOR EQUATION 2 

CONVENTIONAL 

SCS 

TECHNIQUE 






SUMMARY OF ERTS-MSS 
ERTS 

MEANS 9 . 
ORBIT 1 

standard 

NUMBER 1 
GROUP 

DEVI 
274 i 
I 

ATI ON .RUNOFF 
ON 4 23 73 

COE FF I C I 

ENTS 


ws 

MSS-4 

MSS“5 

MSS“6 

MSS-7 

SD“4 

SD-5 

SD”6 

SD~7 

COEF1 

COEF2 

COEF3 

205 

31.41 

27.83 

43.28 

23.14 

0.78 

1.87 

1.87 

1.30 

54.4 

6 1 e 

0.039 

207 

31.67 

34.44 

43.76 

21.72 

0.97 

3 o 40 

2.16 

1.32 

75.8 

86. 

0.122 

in 

33.42 

31.90 

47.52 

25.13 

3.58 

7 o 82 

6.53 

4.28 

60.9 

71. 

0.038 

HI 

33.35 

34.25 

47.04 

24.57 

6.37 

8.84 

7*95 

4.71 

58.0 

74. 

0.023 

512 

30.92 

28.68 

40.10 

20.82 

4.24 

3.71 

5.97 

4.24 

67.2 

74. 

0.050 

513 

31.02 

29.02 

38.80 

19.98 

4.30 

3.40 

5.15 

3.57 

65.7 

74. 

0.054 

5141 

31.39 

29.10 

37.54 

19.22 

5 s 07 

2.35 

3.46 

2.23 

61.5 

74. 

0.041 

5146 

32 © 46 

29.04 

38.13 

19.68 

8.62 

2 o 2 5 

3.55 

2.23 

63.8 

73® 

0.068 

522 

31.99 

28.23 

44.42 

23.58 

6*41 

5« 78 

5 « 42 

3.78 

57.1 

73. 

0.031 

612 

3 0 o 2 9 

27.71 

43.57 

23.06 

1.63 

GROUP 

4.12 
I I 

5.00 

3 o 76 

66.7 

74® 

0*057 

WS 

MSS“4 

MSS- 5 

MSS-6 

M5S-7 

SD-4 

SD-5 

SD“6 

SD“ 7 

COEF1 

COEF2 

COEF3 

206 

31.48 

28.00 

42.48 

22.91 

0 c 9 1 

1 ®37 

1.87 

1 o 38 

5 3 o 6 

61 . 

0*034 

206 

31.26 

30.84 

39.05 

19.89 

Qs73 

3.02 

3 e 36 

2 o 21 

77*4 

8 3 e 

0.147 

121 

34 o 96 

37.86 

50 o 41 

26.08 

7 a 45 

11.04 

9.07 

5 a 41 

58.6 

78. 

0.023 

311 

31 086 

28.19 

49.83 

27.39 

7 o 6 9 

7.61 

6 a 34 

4 o 93 

69 o 6 

77 s 

0.078 

511 

30 o 6 3 

27.40 

43.66 

23.26 

5 © 84 

5 s 2 1 

6.19 

4 e 73 

69 o 4 

75 o 

0.082 

5142 

31 o 1 1 

28.90 

36.51 

18.51 

1 o09 

2 a 71 

3 o 1 0 

2.01 

59 o4 

7 6 o 

0.027 

5143 

3 1 o 1 6 

29.31 

36 s 54 

1 8 o 49 

1 a 10 

2 o 1 1 

3 .01 

1 o 92 

56 o 3 

68 o 

0.021 

5144 

32 o 30 

2 9 & 2 2 

37 o 60 

1 9 o 26 

8 o 5 9 

2 o 3 9 

3 o 43 

2 o 20 

6 2 o 8 

76 o 

0 e 066 

611 

30 a 0 1 

27 o 08 

40 o 2 8 

21 s02 

1 o 6 2 

3 o 96 

5 o 6 3 

3 o 99 

70 o 2 

77 o 

0 o 06 5 

621 

WS 

30 o 66 

OClDa 

30 o 36 41 o 1 3 2 1 o 00 

watershed number 

1 o 56 

4 o 40 

5 o 3 1 

3 o 52 

6 7 o 4 

7 7 o 

0 o 05 7 

COEFl 
COEF2 — 
COEF3 =<=■ 

SCS CURVE NUMBER CALCULATED 
SCS CURVE NUMBER DERIVED BY 
COEFICIENT FOR EQUATION 2 

FROM MEASURED 
CONVENTIONAL 

i RAINFALL AND 
SCS TECHNIQUE 

RUNOFF 




SUMMARY OF ERTS-MSS MEANS eSTANDARD DEV I AT I ON * RUNOFF COEFFICIENTS 
ERTS ORBIT NUMBER 1400 ON 8 2 7 73 
GROUP I 


WS 

MSS-4 

MS 5“ 5 

MSS-6 

MSS-7 

SD“4 

SD*“5 

SD-6 

SD-7 

COE F 1 

CQEF2 

C0EF3 

205 

37*40 

31*60 

46.30 

24.06 

0.77 

2*80 

1*39 

0.76 

54.4 

61. 

0*039 

207 

38 e 74 

36*74 

50.58 

24.63 

1.19 

3.63 

2*32 

0.76 

75.8 

86.' 

0.122 

111 

39 o 56 

37.20 

48.52 

24.07 

2.08 

5.93 

4.41 

2.86 

60.9 

71. 

0.038 

141 

37*09 

32.41 

47.61 

24.20 

2*08 

5.48 

4 o 8 1 

3.01 

58.0 

74. 

0.023 

512 

37*10 

33.15 

45.82 

22.86 

1.62 

6.60 

3.28 

2.07 

67.2 

74. 

0.050 

513 

36 « 92 

33.00 

45.27 

22.55 

1.34 

6.53 

2.85 

1.75 

65.7 

74 * 

0.054 

5141 

37*16 

32.44 

45.23 

22*61 

1.41 

6*68 

2.26 

1.34 

6 1 * 5 

74 o 

0.041 

5146 

37*02 

32.11 

45.54 

23.16 

1.86 

9.22 

2.23 

1.32 

63.8 

73. 

0.06S 

522 

3 8*61 

33*23 

49.84 

25.31 

1*35 

7.10 

3.67 

2.57 

57.1 

73. 

0.031 

612 

0,00 

0 s OQ 

0.00 

0.00 

C.OQ 

0.00 

0.00 

0.00 

66.7 

7 4 « 

0.057 






GROUP 

I I 






WS 

MSS-4 

MSS“5 

MSS-6 

MSS-7 

SD-4 

SD“5 

SD“6 

SC °7 

CQEF1 

C0EF2 

C0EF3 

206 

37*21 

31.75 

46.64 

24.11 

1*13 

2.05 

1.66 

Og 86 

53.6 

61. 

0.034 

208 

37.09 

33.04 

46 * 64 

22*95 

0.92 

1.59 

1.62 

1.17 

77.4 

83# 

0.14? 

121 

37*69 

34.29 

49.71 

25*30 

2.57 

5*68 

7.22 

4.83 

5 8.6 

78. 

0.023 

311 

40.08 

39.44 

46.80 

22.69 

2 e 24 

7*00 

3.95 

2.67 

6? o 6 

77. 

0 e C 7 0 

511 

37*93 

36.34 

47.06 

23.30 

1.64 

7.40 

4.03 

2.68 

69*4 

75. 

0*082 

5142 

37.10 

31.32 

44.64 

22.34 

1.50 

2.60 

2.14 

1*14 

5 9.4 

76* 

0.027 

5143 

36.94 

31.64 

44,40 

22.22 

1 o 1 8 

2*95 

2 o 1 0 

1.25 

5 6 o 3 

6 8 o 

0*02 1 

5 144 

37*01 

32.36 

45 .38 

22 * 86 

1 o 5 5 

9*11 

2 a 19 

1.33 

62 o 3 

76 « 

0*066 

611 

0.00 

OoOO 

0.00 

OoOC 

OoOO 

0.00 

OoOO 

0*00 

70*2 

77 « 

0.065 

621 

38 © 39 

35 o 0 1 

46 o 1 3 

22*73 

2 © 17 

6.90 

4.27 

2.64 

6 7 a 4 

77 o 

0*057 


WS WATERSHED NUMBER 

COEF1 SCS CURVE NUMBER CALCULATED FROM MEASURED RAINFALL AND RUNOFF 
C0EF2 -*= SCS CURVE NUMBER DERIVED BY CONVENTIONAL SCS TECHNIQUE 
C0EF3 COE F I C I ENT FOR EQUATION 2 


> 


& 



SUMMARY OF ERTS-MSS MEANS ^STANDARD DEV I AT ION & RUNOFF COEFFICIENTS 


ERT 5 ORBIT NUMBER 1508 ON 12 13 73 
GROUP I 


WS 

MSS— 4 

MSS— 5 

MSS -6 

MSS-7 

SD-4 

SD“ 5 

SD -6 

SD“7 

CQEF1 

CCEF2 

COEF3 

205 

19*15 

20.31 

23.15 

11.69 

1 oQl 

1.81 

2.17 

0 o 88 

54 e 4 

6 1 e 

0.039 

207 

19.86 

22.36 

24.82 

12.41 

0.99 

2.15 

2 o 10 

1.30 

75.8 

86 . 

0.122 

111 

19*20 

1 8 o 8 6 

22.94 

11.36 

2.07 

3g75 

5.12 

2.95 

60.9 

71c 

0.038 

141 

18.22 

17.41 

20.34 

10.46 

1.84 

3.56 

4.99 

2.88 

58.0 

74 9 

0.023 

512 

17.76 

17.30 

19.88 

10.20 

1.32 

2.64 

4.29 

2 . 46 

6 7 o 2 

74 a 

0.050 

513 

17.65 

17.08 

19.37 

9.92 

1.22 

2.44 

4.02 

2.32 

65.7 

74. 

0.054 

5141 

17.62 

17.11 

13.80 

9.64 

1.15 

2.33 

2.78 

1.49 

61.5 

74. 

0.041 

5146 

17.25 

16.38 

IT. 75 

9.13 

1.05 

1.99 

2.64 

1.43 

63.8 

73. 

0.068 

522 

18.85 

18.37 

21.64 

11.17 

2.06 

3.70 

5.06 

2.80 

57.1 

73. 

0.C31 

612 

17.85 

17.24 

19.81 

10.21 

1 « 33 

2.50 

3.79 

2 o 09 

66.7 

74 . 

0.057 






GROUP 

II 






WS 

MSS-4 

MSS-5 

MSS -6 

MSS-7 

SD-4 

SD“5 

SD -6 

SD°*7 

COEF1 

COEF2 

CQEF3 

206 

18.78 

19.04 

21.86 

11.21 

0.74 

1.82 

1.56 

0.88 

53.6 

61. 

0.034 

208 

18.11 

19.07 

20.89 

10.59 

1.12 

2.28 

2.42 

1.31 

77 o 4 

83 a 

C. 147 

121 

18.49 

18.22 

21.79 

11.17 

2.12 

4 . 4 1 

5.50 

3.05 

58.6 

78 & 

C.023 

311 

18.75 

18.23 

23.38 

12*13 

1.35 

2>66 

4.01 

2 « 46 

69*6 

77. 

0.078 

511 

18.75 

18.60 

22.91 

11.83 

1 a 63 

3.27 

4.84 

2.86 

69.4 

75. 

0.082 

51^2 

17.62 

17.42 

19.10 

10*01 

0.97 

2 9 09 

2.29 

1.30 

59.4 

7 6 . 

0.027 

5 143 

17.81 

17.72 

19.49 

10.08 

1.04 

2.27 

2.61 

1.36 

56.3 

68 . 

0.021 

5144 

17.32 

16.56 

18.02 

9.22 

1.03 

2.03 

2 a 5 9 

1.38 

62.8 

76. 

0.066 

611 

18.01 

17.50 

19.82 

10.20 

1.40 

2.50 

4.08 

2.42 

70.2 

77. 

0.065 

621 

18.24 

18.53 

20.71 

10.48 

1.57 

3.27 

4.29 

2.30 

67.4 

7 7 o 

0.057 


WS — — WATERSHED NUMBER 

COEF1 — SCS CURVE NUMBER CALCULATED FROM MEASURED RAINFALL AND RUNOFF 
COEF2 — SCS CURVE NUMBER DERIVED BY CONVENTIONAL SCS TECHNIQUE 
COEF3 — COEFJCIENT FOR EQUATION 2 



MERGE 



INTEGER RSKIP*TAPID 

DIMENSION I DD <36001 * ID <1800 ) * 10 < 1800 1 * I H 1800 1 
DIMENSION I TABU 5) • IDSBFU1) * I OTBL ( 3 ) * I M < 5 ) *IDENT<13) 
EQUIVALENCE < I DD ( 1 ) » I D ( 1 ) ) • ( I DD ( 1 801 ) * 1 0 ( 1 ) ) ♦ ( I M < 5 1 * I DD ( 1 ) ) 
DATA IDENT/-lfOi 4*5*6*7* *ER » * 'TS 1 * * A */ 

DATA ITABL/'E ‘ * * UC 1 * • BR ♦ * 5 * 1 / 

READ< 2*10) RSKIPi NOEFT»NOELT* NRTCY» TAPI D 
10 FORMAT < 1 A I 5 ) 

CALL DS0R< ITABL15) •IDSBF(ll) * I COMP) 

I DENT ( 2 1 =TAP I D 
IM(A)=TAPID 
I M < 2 ) = I DENT ( 7 ) 

I M ( 3 ) *1 DENT ( 8 ) 

I OTBL ( 3 ) =1 
I OTBL { 2 1=0 
I DENT ( 1 0 ) - RSK I P 
I DENT (11) =NOEFT 
I DENT < 1 2 ) =NOELT 
IDENT < 13 ) =NRTCY 
LL'NOEFT+NOELT 
I OTBL ( 1 ) *- ( LL+LL ) 

DO 2 1=1*2 
J * I “ 1 

DO 2 K= 1 * R SK I P 

CALL MGTAP ( J * 2* IO*1800*N01 

IF<NO)2*200*2 

2 CONTINUE 

DO 20 1 = 1 * NR TCY 

CALL MGTAP<0*2*IO*1806*NO) 

IF( NO 13*100*5 

3 NO® 1648 

5 NO=NO-28 
M=NO-NOEFT+l 
CALL MOVE < I0*M*N0*ID*1) 


> 



MERGE CONTINUED 



CALL MG TAP ( 1 *2#IL*18QO*NO) 

! F ( NO 515*100*15 
15 M=N0£FT+1 

CALL MQVE( I L » 1 #NOELT • I D »M ) 

CALL WDDSK ( I»ID(LL) »I0TBL(3 } * IDSBF( 11) * I COMP ) 
20 CONTINUE 

100 CALL MGTAPC0*6tILt9*N0) 

CALL MGT AP < 1 f 6 * I L * 9 *N0 ) 

DO 2 5 I*ltRSK2P 

CALL MGTAP(0»2t ILtlSOOfNO) 

IF(NC)25*200»25 
25 CONTINUE 

27 CALL MGTAP{l*3fIDENT*13*NO) 

IF(NO)27#200i28 

28 CONTINUE 

DO 50 1*1 » NR TCY 

CALL MGTAP(0t2*IL*l80O#NO) 

IF { NO ) 2 9 • 200 i 30 

29 NO* 1648 

30 NO * NO-2 8 

CALL RDDSK ( I * I DD ( LL ) » IOTBL ( 3 ) t IDSBF( 11) » I COMP) 
K=NOEFT + NQ+1 

CALL MOVE ( IDD*MfLL*IDD»<) 

CALL MOVE ( I L * 1 *N0 * I DD *M ) 

K = LL + NO+4 
I M ( 1 ) = I 

35 CALL MGT AP ( 1 * 3 * I M >K*NO) 

I F ( NO ) 3 5 # 2 00 * 50 
50 CONTINUE 

CALL MGTAP{ 1 »0tIDD#9*NO) 

200 CALL MGTAP{0#6#IDD#9iNO) 

CALL MGTAP ( 1 *6# IDD»9*N0) 

CALL EXIT 
END 



OF POOR QUALETC 


MGTAP 


******** ************ a****************************** 


* IBM 1800 MAGNETIC TAPE I/O ROUTINE * 

* * 

* THE FOLLOWING ROUTINE WILL ALLOW THE USER TO * 

* READ BINARY TAPES IN FORTRAN OR ASSEMBLER. * 

* THIS ROUTINE HAS BEEN SUCCESSFULLY USED TO READ * 

* TAPES CREATED. ON IBM 360»S* DEC 10»S* AND UN I VAC* 

* 1100 SERIES COMPUTERS. * 

* * 

* AUTHOR — * 

* ROBERT J TORLINE * 

* USDA-ARS * 

* P 0 BOX 267 # 

* WESLACO# TX 78596 * 

* 512-968-5533 EXT 53 * 


* COMMON USER NO 3080 * 

* ' * 

* MAJOR REVISION VER 4 DATE OF LAST CHANGE * 

* MINOR REVISION MOD 1 FEB 22* 1974 * 

****************** ******************#.** ## * ## **. w . # 3 ## 

ENT MGTAP ( I UNI T » I FUNC # I D t MAX #NO ) 

MGTAP DC *-* 

ST X 1 . XR 1 + 1 

LDX II MGTAP 

LD II 1 OPERATION CODE 

STO MGTAP 

SLA 12 

OR 11 0 DRIVE NUMBER 

STO CONTR CONTROL WORD 

LD II 3 VALUE OF RECORD LENGTH OR 

STO RL MAX DIMENSION OF I/O ARRAY 

LD 12 ADDRS OF END I/O TABLE 
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MGTAP CONTINUED 



STO 


NAME+1 


S 


RL 


STO 


IOADD 


LD 

1 

A 


STO 


RETRN+1 


STS 

I 

IOADD ♦ / AO 


LD 

I 

IOADD 


STO 


SAVE 


LD 


RL 


STO 

I 

IOADD 


MDX 

1 

5 


S.TX 

1 

BACK* 1 


LD 


MGTAP 


A 


M3 


BN 


SKIP 


A 


Ml 


BP 


SKIP 


LD 


RL 


BSI 


FLIP 


LD 


M3 

SKIP 

STO 


MGTAP 


CALL 


P 2 AO 1 

CONTR 

DC 


*-* 

IOADD 

DC 


*-* 

RETRN 

STO 

L 



BNP 


TEST 


LD 


MGTAP 


BNN 


XR1 


LD 

I 

RETRN+1 

FLOP 

BS I 


FLIP 

XR1 

LDX 

LI 

*-# 


LD 


SAVE 


STO 

I 

IOADD 

BACK 

BSC 

L 

*-# 

SAVE 

DC 


# — » 


ADDRS BEGINNING I/O TABLE 
ADDRS OF ACTUAL REC LENGTH 
OR COMPLETION CODE 
UNSTORAGE PROTECT BEGIN 
OF I/O TABLE 
AND SAVE CONTENTS 


BRANCH ON READ NEG 

BRANCH ON CONTROL POS 

FLIP ARRAY BEFORE WRITE 
NEG FOR WRITE OR READ 
POS OR ZERO FOR CONTROL 


STO COMPLETION CODE 
IF NOT CORRECT WORD COUNT 
TEST* IF READ OR WRITE 
SET UP TO FLIP ARRAY 


RESTORE SAVED WORD 



MGTAP CONTINUED 


RL 

DC- 


*-* 

M3 

DC 


-3 

Ml 

DC 


-1 

FLIP 

DC 


*-* 


STX 

2 

XR2 + 1 


STX 

3 

XR3 + 1 


STO 


COUN T+l 


LD 


RL 


SRA 


1 


CMP 


COUNT + 1 


MDX 


NAME 


STO 


CO UN T+l 

NAME 

LDX 

LI 

*-# 


LDX 

12 

IOADD 

COUNT 

LDX 

L3 


LOOP 

MDX 

2 

1 


LD 

1 

0 


XCH 




LD 

' 2 

0 


STO 

1 

0 


XCH 




STO 

2 

0 


MDX 

1 

-1 


MDX 

3 

-1 


MDX 


LOOP 

XR2 

LDX 

L2 


XR3 

LDX 

L3 

*-* 


BSC 

I 

FLIP 

TEST 

LD 


MGTAP 


S 


M3 


BN2 


XR1 


LD 


RL 


MDX 


FLOP 


END 




ROUTINE TO FLIP ARRAY 


FOR 360 CQMPATABILITY 
DATA ARRAYS MUST BE 
FLIPPED AFTER READS AND 
BEFORE AND AFTER WRITES* 
NOTE — 

- IF AN ERROR OCCURS ON - 

- A WRITE COMMAND THE 

- DATA ARRAY WILL RETURN - 

- TO THE CALLING PROGRAM - 

- AS IT WAS BEFORE THE 

“ CALL* RJT - 


ON ERROR TEST FOR WRITE 
COMMAND {READ WILL BE A 
-I OR -2)* IF WRITE SET UP 
TO FLIP THE ARRAY AND 
RETURN IT AS BEFORE* 
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P2401 



ENT 


P2401 



P 240 1 

DC 



CALL 

P2401 


STX 

1 

XR1 + 1 

DC 

/HEX CONTROL 


LDX 

11 

P 240 1 

DC ' 

ADDRS 


LD 

1 

0 




TTO 


CONTL 




LD 

1 

1 

ACCUM 

+ WORD COUNT 


S TO 


ADDRS 


0 EOF 


MDX 

1 

2 


-1 WLR 


STX 

1 

BACK* 1 


-2 R/W ERROR 

CALL 

CALL 


MAGT 

2401 

-3 END OF TAPE 


DC 


LIST 


-4 FILE PROTECT 


MDX 


* + 2 


-5 REQUEST ABOR 

8USY 

LD 


* 




DC 


/ 3 04 5 




LD 


LIST 




SSC 

L 

BUSY 9 Z 




LD 


LIST+6 




A 


Ml 




BZ 


ZERO 




A 


Ml 




A 


Ml 




BNP 


WAIT 




A 


Ml 




BZ 


EOF 




A 


Ml 




BZ 


WLR 




EOR 


Ml 

CHANGE SIGN OF ERROR 


MDX 


ZERO+2 



WLR 

LD 


L I 5T+2 




BNN 


ERROR 




A 

I 

ADDRS 




S 


Ml 




MDX 


ZERO+2 
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P2AQ1 CONTINUED 





ERROR 

LD 



MDX 


ZERO 

LD 

I 

XR1 

LDX 

LI 

BACK 

BSC 

L 

WAIT 

SLA 



STO 

L 


LD 

WAIT 

MDX 

MDX 

MDX 

L 

SPAUS 

EQU 


EOF 

EQU 


Ml 

DC 


LIST 

DC 



DC 



BSS 



DC 


CONTL 

DC 


ADDRS 

DC 



END 



MI 

ZERO+2 

ADDRS 

*-* 

16 

SPAUS 

CONTL 

SPAUSiO 

WAIT 

CALL 

97 

ZERO+2 

-1 

0 

0 

A 

*-* 

*-# 
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MOVE 



ENT 


MOVE MOVE SUBROUTINE ENTRY POINT 

* 



CALL MOVE! JCARD. Jt JLAST .KCARD *K ) 

* 



THE WORDS JCARD(J) THROUGH 

# 



JCARD(JLAST) ARE MOVED TO KCARD 

* 



STARTING AT KCARD(K). 

MOVE 

DC 


*-* ARGUMENT ADDRESS COMES IN HERE 


STX 

1 

SAVE 1 + 1 SAVE I R 1 


LDX 

11 

MOVE PUT ARGUMENT ADDRESS IN IR1 


LD 

1 

0 GET JCARD ADDRESS 


S 

11 

2 SUBTRACT JLAST VALUE 


STO 


LD1+1 PLACE A DDR OF JCARD ( JLAST ) IN 

# 



PICKUP OF MOVE 


LD 

11 

2 GET JLAST VALUE 

ONE 

S 

11 

1 SUBTRACT J VALUE 


BSC 


+2 CHECK FIELD WIDTH 


SR A 


16 NEGATIVE - MAKE IT ZERO 


STO 


LDX + 1 STORE FIELD WIDTH IN LD < 


LD 

1 

3 GET KCARD ADDRESS 


S 

11 

4 SUBTRACT K VALUE 


s 


LDX+ 1 SUBTRACT FIELD WIDTH 


STO 


STO+1 PLACE ADDR OF KCARD ( KLAST ) IN 

* 



STORE OF MOVE 


MDX 

L 

LDX+ltl ADD ONE TO FIELD WIDTH 




MAKING TT TRUE 


MDX 

1 

5 MOVE OVER FIVE ARGUMENTS 


STX 

1 

DONE 1 + 1 CREATE RETURN ADDRESS 

* 



JNOW-J 

* 



KNOW =K +J NOW- J 

LDX 

LDX 

LI 

*-* LOAD I R 1 WITH FIELD WIDTH 

# 



KCARD ( KNOW )= JCARD < JNOW) 

LD1 

LD 

LI 

*-* PICKUP JCARD ( JNOW ) 

STO 

STC 

LI 

*-* STORE IT IN KCARD ( KNOW ) 

* 



SEE IF JNOW IS LESS THAN JLAST. 

* 



IF YES > JNOW = JNOW+ 1 AND MOVE 



MOVE CONTINUED 


# 


* 

SAVE1 

DONE1 


MDX 

MDX 

1 

-1 

LD1 

NEXT CHARACTER • IF 
DECREMENT THE FIELD 
NOT DONE - GET NEXT 

CVIT 

NO • EX I T • • • • 

WIDTH 

WORD 

LDX 

BSC 

END 

LI 

L 

*-* 

*-* 

DONE - RESTORE IR1 
RETURN TO CALLING PROGRAM 
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OKLAH 


c ##***##***#*# ***************** ************************************ 
C 

C OKLAH 

C PROGRAMMED BY M* GAUTREAUX 

C USDA • WESLACO » TEXAS 78596 

£*******#*************##************************************************ 
C OKLAH IS A PROGRAM DESIGNED TO SELECT IRREGULAR SHAPE 

C AREAS BY A SERIES OF TRAPEZOIDS USED TO APPROXIMATE THE AREA 

£*#*#** a***************************************************** *********** 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I IN 
ISAVE 

IDSPY 

I OUT 

QUAD 

ID 

IPT 

DENT 

JTAPE 

ITAPE 

IFILE 

NODAT 

NOSET 

NRCD 


INPUT VECTOR 

VECTOR USED TO SAVE WANTED DATA 
(DIMENSION MAY VARY DEPENDING ON AREA SIZE) 

VECTOR USED IN DISPLAYING AREA ON DICOMED MODEL 36 
DISPLAY UNIT 

VECTOR USED TO CONTAIN TWO DATA PTs* CONTAIN IN ONE 


ERTS WORD 

VECTOR (TWO DIMENSIONAL) USED TO HOLD THE ID OF AREA 
AND OTHER PARAMETERS FOR THE SELECTED AREA 
ID VECTOR (512) 

VECTOR CONTAINS THE FOUR PTS. THAT DEFINE THE AREA 
( ONE TRAPEZOID AREA) 

VECTOR USED TO CONVERT ID TO REAL NUMBERS 
SECONDARY TAPE ID NUMBER ON TAPE (READ FROM TAPE) 
SECONDARY TAPE ID NUMBER ON TAPE (READ FROM CARD) 

CURRENT NUMBER OF FILES ON SECONDARY JAPE crT cr> 

NUMBER OF PASSES (REWINDS) NEEDED TO GET ALL SELECTED 


AREAS 

NUMBER OF AREAS TO BE SELECTED ON CURRENT PASS 

NUMBER OF RECORD THAT YOU WANT TO START WITH 

WILL SKIP NRCD-1 RECORDS ON ERTS TAPE 

USED MAINLY WITH ERTSA - BY R. TQRLlNE - TO TRANSFER 


& 



^»or» onnnnor>nnnnr>nnnnnnr>nn 


OKLAH CONTINUED 


THE Y STARTING POSITION 

IWORD - NUMBER OF WORDS (ONE WORD ■ FOUR CONSECUTIVE ERTS 

WORDS WHERE THE FIRST ERTS WORD WAS FOR CHANNEL ONE) 
USED WITH ERTSA - BY R. TGRLINE - TO TRANSFER 
THE X STARTING POSITION 
IFLNO - FLIGHT NUMBER OF CURRENT ERTS TAPE 
I H20 - WATERSHED NUMBER 

IMNTH - DATE OF PASS BY SAT. 

I YEAR “ — — - 

IHSTA - BOUNDARY PTS. OF THE LEFT - HAND SIDE OF SELECTED 
AREA (IHSTA * THE LOWEST X COORD.- X) 

IHEND - BOUNDARY PT. OF THE RIGHT - HAND SIDE OF SELECTED 
AREA (IHEND * THE HIGHEST X COORD. +1) 

I X * I Y - CORRECTION FACTORS FOR A OVERALL SHIFT 
M - BLOW-UP FACTOR (FOR DISPLAY PURPOSES) 

I H * I V - STARTING COORD. ON DISPLAY 
ICHAN - CHANNEL (1-4) SEND TO DISPLAY 




INTEGER START 

DIMENSION I INC 1800) *ISAVE<3000) #IDSPY(1GQ0) #IOUT(3600) 
DIMENSION QUAD (60 *12) * I D ( 1 0 ) * I PT ( 8 ) * E ( 4 ) * PT ( 8 ) *DENTl 10) 
OEFINE FILE 1 ( 200 * 8 *U . KAPPA ) 

KAPPA*1 

READS IN JTAPE AND IFLE FROM SECONDARY TAPE 

READ ( 5 ) JTAPE * I FI LE 
WRITE(3*6) JTAPE*IFILE 
6 FORMAT! 2110) 

C 

RE AD (2*767) I TAPE t NO DAT * NO SET *NRCD* IWORD 
767 FORMAT( 515 ) 
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OKLAH CONTINUED 


C 

KKK-0 

MAX-1800 

C 

C CHECKS FOR CORRECT SECONDARY TAPE 

C 

IF (JTAPE-I TAPE) 600*601 #600 

601 CALL W1322(5*2) 

I F ( IFILE*- 1)610*620* 620 
620 LF I LE* I F I LE + 1 
C 

C SKIPS FILES ON SECONDARY TAPE IF NECESSARY 

C 

CALL W1322 ( 5 * LF I LE ) 

NFILE-IFILE 
GO TO 602 
610 NFILE-0 
C 

C #*** ###**#****#****#****#*#*#***##*#** ** + ***#*##****■*■&***##■##■#**##* # 

C *»»*#***#*** ***■-**# ********************* *#■***+*■***##-*■»-■«■#****•** ********** 

C ***************************************************************** ****g 

c 

C CONTROLS NUMBER OF PASSES NEEDED BY NODAT 

C 

C#**##*###*#****#******** **************** ******************************* 

c 

602 DO 630 I2*1*N0DAT 
CALL W1 322 ( 4 * 2 ) 

C 

C READS 2ND PARAMETER CARD - NOSET 

C 

RE AD 12*768) NOSET 
768 FORMAT ( 15) 

C 

C#***####***##*#******* ********* *************** ************ ************* 
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oklah continued 



c 

C CONTROLS NUMBER OF AREAS IN CURRENT PASS BY NOSET 

C 

C ##**## ###*#*#**# »#****««#**#*#*#«**#**«# #***************#*#****#-##****.* 

DO 500 I B IG=1 *NQSET 
C 

I RECD*Q 
C 

C READS 3TH PARAMETER CARD - I NFLNO • IH2Q * I DAY t IMNTH * I YE AR 

C 

READ < 2*237 } IFLNO* IH20* I DAY ♦ IMNTH* I YEAR 

237 FORMAT (515) 

C 

I I N { 1 ) « I FLNO 
I IN ( 2 ) * I H20 
1 1 N ( 3 ) * I DAY 
IIN(4)=IMNTH 
I IN ( 5 } = I YEAR 
DO 238 I A*6 * 9 

238 I X N ( I A ) *0 
C 

C WRITES 3TH PARAMETER CARD ON 1ST RECORD ON SECONDARY TAPE 

C FOR THE NEW FILE 

C 

CALL MGTAP{ 1*3*1 IN*9*NO) 

NRCD*NRCD-1 
DO 222 JN* 1 * NRCD 
222 CALL MGTAP{0*2*IIN*MAX*NO ) 

C 

C# **###*##■*## a*************************** *##« **« **•*■* ***■*•#■***# **#*##*•**** 

c 

C READS 4TH PARAMETER CARD - 

C 

READ(2*100)NQUAD*IHSTA*IHEND*IX*IY*M*IH*IV*ICHAN 
WRITE(3*100)NQUAD*IHSTA*IHEND*IX*IY»M*IH* I V* I CHAN 
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OKLAH CONTINUED 



100 FORMAT (915) 


C 

c 

£4 444444444444444 4 #*##*#*******##*#**#*■**•##*•#■■#■■*#■***####**■*■& #*•*■**# <HHH* * fc * 


c 

c CONTROLS NUMBER OF TRAPEZOIDS IN AREA BY NQUAD 

£4444444444444444444*4*tt#**tt***#********tt******************************* 

C 

DO 200 KK=1 • NOUAD 


C 

C READS IN NQUAD GROUPS OF FOUR RTS • THAT DETERMINE 

C EACH TRAPEZOID THAT WILL COMPRISE THE IRREGULAR SHAPED 

C ARE PLUS AN ID NUMBER FOR EACH TRAPEZOID 

C 

READ 12*101} ( ID 00 #K=1*5) #UPT< I)»I«1»8J 

WRITE(3»101)UD<K)*K*1*5)*UPT(1>#I*1*8) 

101 FORMAT (512*815) 

C 

C SHIFTING OF COORD**S DONE HERE IF NECESSARY 

C 

I HSTA= I HST A+ 1 X 

ihend-ihend+ix 
DO 10 J*1 *8*2 
JK* J+l 

IPT( J)-IPT( JJ+IX 
10 IPT( JK)*IPT( JDM-IY 
C 

C CHANGING ID NUMBER FROM INTEGER TO REAL 

C 

DO 20 JJ = 1 *8 

PT ( JJ ) "FLOAT ( I PT C J J ) ) 

20 CONTINUE 


C 

C CHECKS FOR 1 •0« ♦ DEMON# 

C 
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OKLAH CONTINUED 



A=PT<5>-PT<3) 

I F ( A ) 21 *25»21 
25 A® A+l 

21 B*PT { 1 ) -PT ( 7 ) 

IF(B)23»22#23 

22 B a B+l 
C 

C CALCULATES THE Y INTERCEPTS AND SLOPES OF LINES THAT MAKE 

C UP THE SIDES OF THE TRAPEZOID 

C 

23 E(1)=IPT<6)-PT(A) )/A 

E ( 2 ) ®-( E ( 1 ) *PT ( 3 ) H-PTI4) 

E(3)=(PT(2)-PTI8) )/B 
E(4)®-<E(3)*PT<1) ) + PT ( 2 ) 

C 

C STORES ABOVE INTO QUAD 

DO AO KL S 6 # 9 
L a KL"5 

40 QUAD ( KK tKL ) a E ( L ) 

C 

C STORES I D 1 S INTO QUAD 

C 

DO 30 K J S 1 • 5 

DENT { KJ ) =FLOAT( ID(KJ) ) 

30 QUAD! KK • KJ ) ®DENT ( KJ > 

C 

C STORES PTS« 2 AND 6 INTO QUAD (PTSt 2 AND 6 ARE USED TO 

C COMPLETE THE TRAPEZOID - TOP AND BOTTOM ♦ RESP •- 

C 

QUAD(KK#10)®IPT<2) 

QUAD ( KK 1 11 ) ■ IPT ( 6 ) 

200 CONTINUE 
C 

C#****# *************************************************************** *********** 

c 
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QKLAH CONTINUED 




MAX = 1800 
J a 0 
C 

C CLEARS DISPLAY UNIT OF ALL COMMANDS 

C READYS DISPLAY UNIT FOR START OF INPUT DATA 

C 

CALL EXTIN 
CALL SOI 
C 

£##**##*#***#*#####**#*#**##**»********##******#*****###**###■#**#**#•**■****#***** 
C*##*#**###** #*#*#**#**** #*****##******##********###**#******##**#####.#***.*.**.*** 

C 

DO 300 KK = 1 #NQUAD 
C 

C BEGINS BACK THE NEEDED DATA FACTORS FOR EACH TRAPEZOID 

C PROGRAM WORKS WITH ONLY ONE TRAPEGQ I D AT ANY ONE TIME 

C 

IVSTA=IFIX(QUADUKf 10) ) 

I VEND® IFlX(QUAD(KKtll) > 

C 

C POSITIONS AREA ON DISPLAY WITH ADJUSTMENTS MADE FOR IH*IV* 

C AND BLOW-UP FACTOR 

I HORZ* I H+I HSTA*M 
I VERT* I V+ 1 V$TA*M 
CALL RANDM ( I HQRZ # I VERT > 

C 

C 

C CHECKS POINT IHSTA TO SEE IF FIRST PT# IS IN A ODD OR EVEN 

C POSITION IN INTERWEAVED ERTS DATA 

C 

I IB I T* I SLA ( IHSTA* IS ) 

C 

NUMB*A*( IHEND-IHSTA+1) 

NB*A*< ( IHEND-IHSTA+1) /2 ) 

BNUMBaNB 
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OKLAH CONTINUED 




KBEGN=l+4*( { IHSTA-1 )/2 ) 

C 

11 CALL MGTAP ( 0»2 # I IN*MAX »NO ) 

J*J+1 

C 

C CHECKS TO SEE IF THE Y COQRD.(PT*2) HAS BEEN REACHED 

C 

IF( J-IVSTA) lltl5«15 
C 

15 I F (NO ) 16 • 125 1 18 
C 

C CONSIDERATION FOR A BAD READ CONDITION ON 1ST TAPE 

C 

16 WR I TE ( 3 » 104 ) J 

104 FORMAT { 9 BAD RECORD ’ ♦ 14 ) 

GO TO 11 
C 

C CHECKS CHANNEL 4 FOR STARTING OF GOOD DATA 

C 

18 I CH£K*4 

4 IF ( I IN ( I CHEK ) ) 2 f 1 * 1 

1 START 3 I CHEK 
GO TO 51 

2 I CHEK= I CHEK+4 

IF ( ICHEK-1740 >4*4tl2 

12 WR I TE ( 3 # 1 09 ) 

109 FORMAT ( * ALL NEGATIVE*) 

GO TO 16 
C 

c 

C COORNIATES THE X COMPONENT - FOR ERTSA 

C 

51 IF ( IWORD ) 62 *62*63 
63 START = l + ( I WORD-1 ) *4 
62 KBEGN=5TART+4*( ( IHSTA-1 )/2) 
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OKLAH CONTINUED 


o S 









65 N«NO-START-27 
NB*BNUMB 

IF< I HEND—N >72*72*61 
61 WR ITE ( 3 * 1G5 ) 

105 FORMAT ( » IhEND WENT TOO FARM 
GO TO 300 
C 

C SPLITS ERTS WORD TO GET THE TWO DATA POINTS CONTAINED 

C IN ONE ERTS WORD 

C 

72 CALL HALF ( I IN ( K8EGN ) • IOUT *NB ) 

DO 75 I*1*NB 
I F ( I OUT ( I ) ) 73 * 75 *75 
73 1*1+1 
GO TO 77 
75 CONTINUE 
1 = 1 

77 NB-NB-I+1 
C 

C CHECKS FOR EVEN OR ODD POSITION 

C 

IF ( I1BIT)201 >201*202 

201 JM=2 
NN*2 

GO TO 83 

202 JM= 1 
NN = 1 

C 

83 ICONT = 0 
KL*1 
LJ=4 

Y=FLOAT ( J ) 

C 

C #* ******* **###**# ******************* a*#****##**#*#-*** #** * « 

C HERE THE JTH LINE (Y) AND ITH POSITION (X) ARE > 



n> r\ r\ ooo noonoo 


OKLAH CONTINUED 


CHECKED ON A POINT BY POINT BASES IN THE AREA 
DEFINE BY IHSTA* IHENDoPTo 2 * AND PTe6 

ONLY POINTS LYING IN THE AREA DEFINE BY THE TRAPEGOID 
ARE KEPT 9 OTHER POINTS ARE ZEROED OUT 


DO 155 I=IHSTApIHEND 
XeFLCAT ( I ) 

IF( J°IVEND> 19*299*299 
19 IF(J ra lVSTA}4li>32o32 

32 EG2=QUAD(KKe6)frX*QUAD(KK«7)-Y 
I F ( QUAD ( KK p 6 ) >33e34*34 

33 J F t EQ2 ) 35 0 35 s>41 

34 1 F I EQ2 ) 41 o 3 5 o 35 

35 IF (J^IVEND) 36e41e4l 

36 £Q4=QUAD<KK#8I*X*GUAD<KK«9>~Y 
I F ( QUAD ( KK 9 8 ) )37s38«38 

37 I F ( EQ4 ) 41 e 39 « 39 

38 IF(EQ4)39s39*41 

SAVE GOOD DATA 

39 DO 165 K I =KL o L J 
ISAVE t KI) = IOUT ( JM) 

165 JM-JM42 
JM 3 JM“2 
GO TO 144 

ZEROES OUT UNWANTED DATA 

41 DO 166 Kl-KLolJ 

166 ISAVEUI)*0 
JM-JM+6 

C 

144 lCCNT-ICONT-i-4 
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OKLAH CONTINUED 




KL a KL+4 
L J S L J+4 

GO TO(175* 176) tNN 

175 NN*2 
JM* JM-5 
GO TO 155 

176 NN= 1 
JM* JM+1 

155 CONTINUE 
C 

C WRITES ISAVE VECTOR CONTAINING DATA OF TRAPEQO I D ON TO 2ND TAPE 

C 

CALL MG TAP ( 1 1 3 # ISAVEtNUM9 *KACNB ) 

C 

C 

C BLOWING - UP OF ISAVE VECTOR -IDSPY- FOR DIPLAYING 

C 

NM = 0 

DO 180 KM=ICHAN*NUMBt4 

DO 180 

NM-NM+1 

180 IDSPY(NM)=I$AVE{KM> 

C 

C DISPLAYS IDSPY VECTOR 

C 

DO 190 NK=1«M 
CALL D I SPY I 2 DSPY #NM ) 

19/ CALL EOL 
C 

C CHECKS TO SEE IF REACHED END OF SCREEN 

C 

IF( J-2048)170il70»125 
C 

C COUNTS NUMBER OF RECORDS THAT COMPRISE A FILE 

C 
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OKLAH CONTINUED 


170 I RECD* I RECD+1 
C 

GO TO 11 
C 

299 J*J-1 
C 

C DUE TO THE FACT THAT THE BASE OF ONE TRAPEOOI D IS THE TOP 

C OF THE FOLLOWING ONE A BACK-SPAC E I NG OF ONE RECORD IS NEEDED 

C BACKSPACING ONE RECORD DONE HERE 

C 

CALL MGTAP ( 0 # 7 * I5AVE *NUMB # I EOF ) 

C 

300 CONTINUE 
C 

£********#******#************************* ******** * ************************* **** 
c**#*##*#*# ********************* ********************************** ************** 
C END OF FILE MARK COMMAND ON 2ND TAPE 

CALL MGTAP ( 1 f 8 • I SAVE* NUMB t 1 EOF ) 

NFIL£=NF ILE+1 
C 

C SAVES IMPORTANT FACTS ON TEMP* DISK FILE 

C 

WRITE(1 * KAPPA ) JTA PE • IRECD *NOSET * ID ( 1 ) *ID<2)»ID<3)*IH20 
C 

KKK*KKK+1 

C 

500 CONTINUE 
C 

£#######*##****##***********#*#*****#****************************#******#******* 
C##*#*#*##**###***#*###****##** *********** ************* * ********* ****** ********* 

C 

C REWIND COMMAND FOR 1ST TAPE 

C 

CALL MGTAP { 0 * 5 # I SAVE ♦ NUMB # I EOF > 

C 
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noon on non noon no 


OKLAH CONTINUED 


630 CONTINUE 

C****#****#***#***#**#*#************-***##-*#************#### #■-#■*#***■*■»■* ******* **** 
£***■***#■*■** ******* *************************************************** *********** 
£***************** ******* **************************************** *************** 


WRITES FACTS STORED OUT ON TEMP. DISK FILE ON TO PRINTER 
125 KAPPA= 1 

IF(KKK-1)131»6A0*640 
6A0 CONTINUE 

WR I TE ( 3 # 877 ) 

877 FORMAT < 1H1 ) 

DO 622 LL* 1 .KKK 

READ ( 1 * KAPPA ) JTAPE . IRECD.NOSET • X D < 1 ) *IDII) »ID(3) »IH20 
WRITE ( 3 #603 ) JTAPE* IDC 1 > . ID<2) . IDO ) t IH20 »NFI LE# IRECD 
603 FORMAT ( 10X * 1 TAPE NUMBER ISSlSt* AREA OF IDS5I2*' IS IN FILE * 
1 1 NUMBER 1 * I 6 * 1 WITH*. 16. RECORDS*) 

622 CONTINUE 

COMMAND TO DISPLAY UNIT THAT ALL DATA HAS BEEN SENT 
131 CALL STOP 

REWIND UNLOAD COMMAND TO 1ST TAPE 
CALL MGTAP < 0 . 6 1 1 SAVE ,NUMB .KACNB ) 

NFILE*NFILE+1 

CALL MGTAP ( 1 .6 • ISA VE. NUMB. KACNB ) 

UPDATES FILE COUNT ON 2ND TAPE 
IFILE CHANGED 


WRITE(5 ) JTAPE.NFILE 
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OKLAH CONTINUED 


£********#****##****#*****************#***************************************** 

C 

CALL MGTAP ( 1 *5* 1 SAVE# MAX •NO ) 

CALL EXIT 


600 WRITE(3»103) ITAPE , e , 

103 FORMAT! 21X# ‘WRONG TAPE -SHOULD BE TAPE NUMBER '#15) 

GO TO 131 
END 
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W1322 




* 



ENT 


Wl 3 2 2 (LUNiN) 

W 1 3 2 2 

DC 


#-# WHERE 


STX 

1 

XR1+1 LUN 


LDX 

11 

W 1 3 2 2 N 


LD 

11 

0 


S 


CA 


STO 


LUN 


LD 

11 

1 


MDX 

1 

2 


STX 

1 

BACK+l 


BNP 


STO 


S 


Cl 


3Z 


XR1 RETURN 

STO 

STO 


N 


CALL 


FILE 


DC 


LUN 


DC 


N 

XR1 

LDX 

LI 

*-* 

BACK 

BSC 

L 


CA 

DC 


A 

Cl 

DC 


1 

LUN 

DC 


*-* 

N 

DC 


#-* 


END 


= 4 OR 5 
a + FILE MARK-1 
1 RETURN 
IF LARGER SKIP 
THAT MANY FILE 
MARKS 

= 0 UNLOAD DRIVE 
= - NUMB OF RECORD 
TO SKIP 

IF ONLY 1 TILE 
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FILE 




ENT 
DC 
STX 
LOX 
LD 
A 

STO 
LD 
A 

STO 
LD 
STO 
MDX 
STX 
8 Z 
BP 
EOR 
S 

STO STO 
LDXR1 LDX 
CALL 
DC 
DC 

CALL CALL 
CMND1 DC 
DC 
BSC 
MDX 
BZ 
A 

BZ 

LD 

SKP 

MDX 


FILE (I*N) 
*-* WHERE 
1 XR1 + 1 
II FILE 
II 0 
READ 
CMND1 
II 0 

UNLOD 
CMND2 
II 1 
FILE 
I 2 

I B A C K + 1 
UNLD 
STO 
MASK 
MASK 
LDXR1+1 


I« 0 OR 1 (MAG DRIVE) 

N* + NUMB OF FILE MARKS 
0 UNLOAD 

- NUMB OF RECORDS TO 
SKIP (ABSOLUTE VALUE) 


- IF A FILE MARK IS 

- ENCOUNTERED DURING A 

- RECORD SKIP IT IS 

- COUNTED AS A RECORD. 


CONVERT NEG RECORD 
SKIP TO POS COUNT 


LI *-* 

P 1053 TAKE EAC PRINTER 

0 OFFLINE 

LIST 

P2401 


LIST 

-Z 

CALL 

MODFY EOF INCOUNTERED 

C3 IF END OF TAPE 

UNLD UNLOAD DRIVE 

FILE 


CALL 


SKIP IF RECORD SKIP 



FILE CONTINUED 


MODFY 

MDX 


MDX 

UNL 

CALL 

DC 

DC 

XR1 

LDX 

BACK 

BSC 

UNLD 

CALL 

CMND2 

DC 


DC 


MDX 

READ 

DC 

UNLOD 

DC 

MASK 

DC 

C3 

DC 

LIST 

DC 


BSS 


END 


1 -1 
CALL 
P 1 05 3 
/0100 
LIST 
LI *-* 

L *-* 
P2401 
*-* 

LIST. 

UNL 

/2000 

/6000 

-1 

3 

6 

6 


PUT EAC PRINTER 
ONLINE 
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EXTIN 



r %t 



ENT- EXTIN 

EXTIN DC *-# 

CALL WRITE 

DC CODE 


RETRN BSC I EXTIN 
CODE DC /8000 

END 


SEND 

RETURN 

CODE FOR EXTERNAL INITIAL 
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WRITE 




ENT 


WRITE 

ENTER SUBROUTINE 

WRITE 

DC 


*-# 

STORES VALUE OF RETRN-1 


STX 

1 

SAVE 1 + 1 

SAVE XR1 


LDX 

11 

WRITE 

LOAD ADDRESS OF VALUE 


LD 

11 

0 

TO BE WRITTEN INTO ACCUM. 


MDX 

1 

1 

MODIFY RETURN ADDR BY 1 


STX 

1 

RETRN+1 



STO 


DATA 

AND STORE 

BUSY 

XIO 


DSWO 

TEST DO FOR BUSY STATE 


BOD 


BUSY 

LOOP IF BUSY 


XIO 


WRIT 

SEND OUTPUT TO DISPLAY 

SAVE 1 

LDX 

LI 



RETRN 

BSC 

L 


RETURN 


BSS 

E 

0 


DSWO 

DC 


0 



DC 


/6701 


WRIT 

DC 


AREA 



DC 


/65CO 


AREA 

DC 


/A002 



DC 


126 


DATA 

DC 


*-* 



END 
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G tjr 


ENT SOI 

SOI DC *-* 

CALL WRITE 

DC START 

RETRN BSC I SOI 
START DC /SQQ1 

END 


SOI 


SEND 

CODE FOR START OF INPUT 
RETURN 


A3 6 



RANDM 


ENT RANDM 

RANDM DC *-* 

STX 1 SAVE 1 + 1 

LOX II RANDM 



LD 

11 

0 


BN 


VERTI 


CMP 


MAX 


MDX 


VERTI 


NOP 



HORIZ 

SRT 


8 


STO 


POSIT+1 


SLA 


16 


SLT 


8 


STO 


POSIT 


LD 


HOR 


STO 


COMND 

BUSY 1 

XIO 


DSWO 


BOD 


BUSY 1 


XIO 


WRIT 

VERTI 

LD 

11 

1 


BN 


RETRN 


CMP 


MAX 


MDX 


RETRN 


NOP 




SRT 


8 


STO 


POSIT+1 


SLA 


16 


SLT 


e 


STO 


POSIT 


LD 


VER 


STO 


COMND 

BUSY2 

XIO 


DSWO 


BOD 


BUSY 2 


XIO 


WRIT 


LOAD NOR 1 2 POSITION IN ACC 
IF - CHECK VERT POSITION 
TEST HORIZ WITH MAX=2Q33 
IF GREATER BRANCH TO VERTI 
IF LESS BRANCH TO HORIZ 
SAVE AND 

STORE 2ND POSITIONING WORD 
ZERO ACCUM* 

RECOVER AND 

STORE 1ST POSITIONING WORD 
SEND THE 

HORIZONTAL RANDOM COMMAND 
AND THE FIRST AND SECOND 
WORD OF THE RANDOM 
POSITIONING FORMAT 
LOAD VERTI POSITION IN ACC 
AND PREFORM THE SAME TESTS 
AS FOR HORIZ* BUT RETURN 
IF VERTI VALUE IS NOT RITE 
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RANDM CONTINUED 


RETRN 

MDX 

1 

2 


STX 

1 

BACK-4- 

SAVE 1 

LDX 

LI 


BACK 

BSC 

L 

*-# 

MAX 

DC 


20A7 

HOR 

DC 


/ 8008 

VER 

DC 


/8009 


BSS 

E 

0 

DSWO 

DC 


0 


DC 


/6701 

WRIT 

DC 


AREA 


DC 


/6SC0 

AREA 

DC 


/400A 


DC 


126 

COMND 

DC 


*-* 

POSIT 

DC 


■ #-* 


DC 


*-* 


END 
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STOP 


ENT STOP 

STOP DC *-* 

CALL WRITE 

DC EOL 

BSC I STOP 
EOL DC /8Q0 A 

END 
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MEAN4 


§>§ 

gg 

R'e? 


REAL MEAN { A ) *N(4) 

INTEGER TAPE 
DIMENSION I SAVE ( 4000 ) 

DIMENSION SUMX14) *SUMXX(4) *5D(4) 

DIMENSION NFILE(32) 

DATA Xl*X2/»STAN f * ‘DARD* / 

MAX=4Q0 0 
KONST «3 
IBACK^lOO 

READ (2*100) TAPE* IDRIV*NPPB*NUFIL 
WRITE { 3 *100 ) TAPE* I DR I V * NPPB * NUF I L 
READ (2*100) NFILE 
WRITE(3*100) NFILE 
100 FORMAT ( 1615 ) 

RE AD ( I DR I V ) JTAPE*IFILE 
WRITE( 3*102 3 J JTAPE* I FILE 
1023 FORMAT ( IX * 2 I 5 ) 

IF(TAPE-JTAPE)90*10>90 

10 IM-IDRIV-4 

DO 40 IN=1*NUFIL 
M I = I N+l 

CALL FLSRH ( ISAVE*MAX*NFILE( IN) * NFILE (MI) * KONST * I DR I V > I BACK ) 

11 CALL MGTAP (IM*1 * I SAVE »MAX*NO ) 

WRITE (3* 120) USAVEUB) *IB = 1.5> 

120 FORMAT ( 1H1 * IX * * SCENE 1274 '*I6*‘ DOVER **I6* 

1‘ DATA FOR THE , *I6* SITE 3 * # I 3 * 2X * ♦ MONTH OF 19**12) 

3 I CONT »0 

DO 400 1-1*4 
MEAN ( I ) 

NII)=0. 

5UMX ( I ) *0. 

SUMXX ( I )=0* 

400 CONTINUE 

26 CALL MGTAPi I M • 1 • I SAVE * MAX * NO ) 
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AL PAGE IS 


MEAN4 CONTINUED 



I F ( NO ) 90 * 27 * 25 
25 L= 1 
LL = 2 

DO 200 J=1*N0*4 

IF( I SAVE ( L ) + I SAVE (LL))50*13*50 
50 JK«J 

DO 12 11*1*4 
6 I CONT * I CQNT + 1 
N(ID s N( II )+l* 

X* FLOAT ( I SAVE ( JK ) > 

SUMX ( II ) =SUMX ( 1 1 ) +X 
SUMXX( I I)=SUMXX( I I >+x*x 

12 JIC= JK+1 

13 L*L+4 
LL=LL+4 

200 CONTINUE 
GO TO 26 

27 CALL MGTAP ( I M * 7 * I SAVE *MAX *NO ) 

WRITE(3*701)X1*X2 
WR I TE ( 3 *702 ) 

701 FORMAT { /*34X*2A4*6X* ' NUM « OF PTS* 1 ) 

702 FORMAT ( 2 5X * 1 MEAN * * 5X* ' DEVIATION' * 5X * ' CONSIDERED 1 ) 
DO 201 I J* 1 * 4 

MEAN ( I J ) =SUMX ( I J ) /NCI J) 

SD( 1 J)=SUMXX( I J } “ ( SUMX ( IJ )*SUMX( IJ) )/N( IJ) 
SDUJ}*SQRT(SD( IJ)/(N{ I J ) -1 • ) ) 

KN* I J 

WRITE(3*703)KN*MEAN(KN> *SD(KN) *N(KN) 

201 CONTINUE , 

703 FORMAT ( 10X* ‘CHANNEL 1 *13*2X*F7*2*3X*F8*3*6X*F9.0) 

40 CONTINUE 
GO TO 80 

90 WR I TE 1 3 * 1 06 ) TAPE 

106 FORMAT* 10X* ‘WRONG TAPE* SHOULD BE TAPE NUMBER' *13) 
80 CALL MGTAP* I M *6 * I SAVE *MAX * NO ) 

END 
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FLSRH 


SUBROUTINE FLSRH ( 1 SAVE*MAX *NF I LE * KF I LE * KONST * I DR I V * I B ACK ) 
DIMENSION I SAVE ( 1 ) 

I M® I DR I V—4 

CALL MG TAP ( I M*7 * I SAVE tMAX #NO ) 

ICHEK=NFlLE-KFlLE 
IF( ICHEK-NFILE)8*10*8 

8 I F ( ICHEK ) 1 * 1 * 2 
ICHEK=KFILE-NFILE+1 
CALL W1 3 22 < IDRiVtICHEK) 

GO TO 10 

2 IF( ICHEK-K0NST)3#3*4 

3 IF( IBACK-NFILE>9*9*4 

9 DO 6 1=1*2 

6 CALL MGTAP { I M #7 « I SAVE * MAX # NO 5 
CALL MGTAP { IM»2#I SAVE * MAX *NO ) 

IF (NO) 3 *7*3 

7 I CHEK = I CHEK- 1 

I F ( ICHEK } 10*3*3 

4 CALL MGTAP(IM*5»ISAVE.MAX*N0) 

CALL W1 322 ( I DR I V * 2 ) 

IF(KFILE-1)10*10*5 

5 CALL W1322( IDRIV#KFILE) 

10 RETURN 

END 



DESPY 



C O^TH^DICOME^D^SPLAY WILL display DIGITAL Data IN secondary format 

INTEGER COR 

i° I DSM 7 00 ) IVECT(6000 ) * i DSAV(700) * I D ( 26 ) » I DD ( 2 56 > *MINMX<2) 

DATA M1NMX/* S'Z 1 / 

IMAX=6000 
C0R*2 
L INE*3 

C NCHAN IS THE CHANNEL THAT IS TO BE DISPLAYED 
C MG IN T IS THE DRIVE ON WHICH THE INPUT TAPE IS LOCATED 

C NUMCH IS THE NUMBER OF CHANNELS ON THE DATA TAPE TO BE DISPLAYED 

C IMAGE WILL RESULT ’ 

C IRCD1 IS THE FIRST RECORD TO BE DISPLAYED 
C IRCD IS THE NUMBER OF RECORDS TO BE DISPLAYED 
C INCRR - INCREMENT RECORDS 

C KSMP1 IS THE FIRST PIXEL ON THE SCAN LINE TO BE D I r PLAYED 
C KSMP IS THE NUMBER OF PIXELS TO BE DISPLAYED 
C INCRS - INCREMENT SAMPLES WITHIN RECORDS 
C IVAL1 - MINIMUM DIGITAL VALUE EXPECTED 

C IVAL2 - MAXIMUM DIGITAL VALUE EXPECTED 

C INCRV - INCREMENT OF DIGITAL VALUES 

READ <2* 1000) NCHAN »MGIN •NUMCH# I RCDlf IRCD »INCDR«KSMPl»ICSMPt INCRS* 

1IVAL1*IVAL2* INCRV 
READ (2*103) ID 
103 FORMAT ( 26A 1 ) 

DO 5 1*1*256 
5 IDD{ I)»MlNMXm 
DO 3 1=1*1 VAL 1 

3 IDDU )=MINMX(1) 

DO 4 I* I VAL2 *256 

4 IDD ( I ) a MINMX ( 2 ) 

1000 FORMAT (1615) 

C INITIALIZE SCREEN 
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DESPY CONTINUED 


J = 0 

K«IVAL1-1 

DO 1 1= I VAL1 * IVAL2 * INCRV 
JaJ+1 

DO 1 L= 1 1 INCRV 

1 IDD(K)sIDU) 

DO 2 I s I VAL1 * IVAL2 t INCRV 
It**I + INCRV-l 

2 WR I TE ( 3 » 101 } I » 2 1 1 2 DD ( I) 

101 FORMAT { 1X»2I8»2X#A1) 

WRlTE(3tl02) 

102 FORMAT(lHl) 

CALL MG TAP (MG IN * 1 * IVECT * IMAX *NQ ) 

C CHECK FOR AN ERROR 
I F ( NO ) 20 * 20 * 40 
20 WRITE(LINEf2000)MGlN>N0 

2000 FORMAT* IX* 1 BAD READ ON DRIVE 1 * I 3 #2X * 1 ERRCODE * ♦ XX * 1 5 * / / ) 

CALL STOP 
CALL RWIND(MGIN) 

CALL EXIT 
AO CONTINUE 

DO 50 I *1 #NUMCH 
I A* 1+2 

IF ( IVECT ( I A5-NCHAN) 50*70*50 
C PAGE I CONTINUE PROGRAM DESPY 
50 CONTINUE 

WR I TE ( L I NE *20 10 ) 

2010 FORMAT( IX* ’COULD NOT FIND CHANNEL POSITION AS SPECIFIED ON TAPE*'C 
#H£CK AND MAKE SURE THAT YOU HAVE RIGHT TAPE-CALL EXIT 1 *//) 

CALL EXIT 

C THIS IS THE CHANNEL POSITION OF THE PIXEL 
70 I POS I a 1 + 24 





DESPV CONTINUED 


KSMP2 =KSMP 1+KSMP-l 
85 CALL MGTAP(MGIN#1#IVECT # I MAX * NO ) 

IF (NO) 20 #20*90 

C THIS IS IF WE GOT A GOOD READ 
C SKIP TAPE RECORDS UNTIL WE GET TO IRCDI 
90 IF( IVECT< l)-IRCDl>85#110#85 
110 IF<lVECT(2))120fl30*120 
120 IF{IVECT(2)-1>125#130»125 
125 WRITE (LINE #2030 ) 

2030 FORMAT (IK# 1 TAPE RECORDS ARE NOT IN CORRECT ORDER-CALL EXIT 1 *//) 
CALL EXIT 

130 I THIS* I VECT ( 1 ) 

IVCTE»IRCD1+IRCD 
MM*MM+ 1 
ID ISE*0 

142 I D I SS* I D I SE 

CALL SELET U VECT UPOSI ) • I VECT #NUMCH#NO ) 

CALL SPLIT I I VECT# IDSAV(IDISS) #NO) 

IDISE*IDI SS+NO-2 

145 CALL MGTAP(MGJN#1#IVECT * I MAX *NO ) 

IF(NO)20*20#150 
150 IF(IVECTU) )160#145#160 
160 IF( I VECT ( D-ITHIS) 170 #142 #170 
170 IF(MM-MMM) 175*180*175 
180 MMM=MMM+ I NCDR 

DO 155 KK*KSMP1#KSMP2»INCRS 
KKK=701-KK 
JJ* IDSAV ( KKK ) +1 
155 XDSAUK) = IDD< JJ) 

WRITE <3 *100) MM* { IDS A UK) #KK*KSMP 1 *KSMP2 • I NCRS ) 

175 IF( I VECT ( D-IVCTE) 130*130*260 
100 FORMAT ( lX*I5#2Xfl24Al) 

260 CALL SK I P C 3 » 1 ) 

CALL EXIT 
END 
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SELET 



***** **** *$ CALL SELET < ID * I OUT * NV ♦NCONT ) 



ENT 


SELET 

$WK4 

EOU 


54 

SELET 

DC 


*-* 


CALL 


QZSAV 


MDM 


$ WK4 * +4 


LDX 

11 

SELET 


LD 

1 

0 


STO 


LDXR1+1 


LD 

1 

1 


STO 


LDXR2+1 


LD 

11 

2 


M 


(CM INI 


SLT 


16 


STO 


NV1 + 1 


STO 


NV2 + 1 


LD 

1 

3 


STO 


ADRS+1 


LD 

11 

3 


STO 


LDXR3+1 

LDXR1 

LDX 

LI 

*-* 

LDXR2 

LDX 

L2 

*-* 

LDXR3 

LDX 

L3 

*-* 


LD 


XO 


STO 

I 

ADRS+1 

LOOP 

LD 

1 

0 


STO 

2 

0 

ADRS 

MDX 

L 

*-* * + l 

NV1 

MDX 

LI 

*-* 


MDX 

2 

-1 

NV2 

MDX 

L3 

*-* 


MDX 


LOOP 


CALL 


QZEXT 

KMIN1 

DC 


-1 


END 
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D r SCM 




SUBROUTINE P I SCM 

C MULTIPLE GROUP DISCRIMINANT ANALYSIS. A COOLE Y-LQHNES PROGRAM. 

C THIS PROGRAM COMPUTES DISCRIMINANT FUNCTIONS# THEIR CANONICAL 

C CORRELATIONS WITH GROUP MEMBERSHIP DUMMY VARIATES# F-RATIOS FOR 

C THESE# AND CENTROIDS OF GROUPS IN THE STANDARDIZED DISCRIMINANT 

C FUNCTIONS SPACE. COEFFICIENTS FOR COMPUTING STANDARDIZED 
C DISCRIMINANT FUNCTIONS SCORES FROM DEVIATION TEST SCORES ARE 
C PUNCHED OUT. 

C 

C REQUIRED SUBROUTINES ARE DIRNM AND HOW. 

C 

C INPUT 

C 

C 1) FIRST TEN CARDS OF DATA DECK CONTAIN A TEXT DESCRIBING THE 
C JOB# WHICH WILL BE REPRODUCED ON THE OUTPUT. DO NOT USE COLUMN 1. 

C 2) CONTROL CARD ( CARD11 ) 

C COLS 1-2 M=NUMBER OF VARIABLES 

C COLS 3-5 KG=NUMBER OF GROUPS 

C COLS 6-10 N=NUMBER OF SUBJECTS 

C COLS 11-12 KC=NUMBER OF CONTROL VARIABLES PREVIOUSLY 

C PART I ALED OUT BY CQVAR (THIS VALUE WILL BE A 

C ZERO IF INPUT MATRICES COME FROM MANOVA). 

C 3 ) T MATRIX (TOTAL SAMPLE DEVIATION SSCP# AS PUNCHED BY MANOVA 

C OR COVAR). 

C A) W MATRIX (POOLED WITHIN— GROUPS DEVIATION SSCP MATRIX. AS PUNCHE D 
C BY MANOVA OR COVAR). 

C 5) GROUP MEANS (AS PUNCHED BY MANOVA OR COVAR). 

C 6) GRAND MEANS (AS PUNCHED BY MANOVA). 

C 

DIMENSION A{20#20)#B(20#205#C(20#20)*T(20)»U(20}#V(20)# 

* W(20)#X(2Q) #Y120) #Z(20) #TIT(10#2G) #D(2C#20) »E(20#2C) • I WS ( 2 0 } » 

* ISCN(20) #NV( 20) 

COMMON M * KG # N * K.C * T I T #NTITLE* I W S . ISCN#NV 
C 
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DISCRM CONTINUED 



1 WRITE (6*2) 

2 F0RMAT!65H1MULTIPLE GROUP DISCRIMINANT ANALYSIS. A COOLE Y-lOHNE S 

2PROGRAM. ) 

2301 FORMAT! 162!10A4) ) 

DO 3 I « ltNTITLE 

3 WRITE!6#4) (TITU * J) »J=1*20) 

WRITE! 6 * 2302 ) ( NV { K> #K = 1 #M ) 

2302 FORMAT! ///IX # • FOR THIS ANALYSIS THE FOLLOWING DATA WERE USED » 

1//20X**BANDS* •515) 

DO 2303 I = 1 • KG 

2303 WRITE(6*230A) ISCN( I) »IWS( I ) 

2304 FORMAT! lXf *ERTS SCENE'iIS*’ WATERSHED ' 1 1 5 ) 

WRITE!6#2305) 

2305 FORMAT!////) 

4 FORMAT ( 20A4 ) 

9 FORMAT (12* 13* 15# 12) 

5 FORMAT ! 1 OX » 5E14.7) 

DO 6 J-ltM 

READ (3*2301) ! C ! J #K ) #K= J #M ) 

6 CONTINUE 
REWIND 3 

DO 7 J = 1 # M 

READ (4# 2301) !B! J#K) *K=J»M) 

7 CONTINUE 
REWIND 4 

DO 8 J * 1 * M 
DO 8 K=J#M 
C ! K# J ) »C ! J # K ) 

8 B!K#J)«B! J*K) 

DO 15 J«1#M 
DO 15 K a 1 #M 

15 A! J*K)-C! J#K> - B!J»K) 

C A NOW CONTAINS THE A MATRIX ( AMONG-GROUPS DEVIATION SSCP MATRIX). 

C B CONTAINS THE W MATRIX ( WI TH I N-GROUPS DEVIATION SSCP MATRIX). 

IF ! M-KG ) 10*11#11 
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DISCRM CONTINUED 



10 MD=M 

GO TO 12 

11 M0=KG-1 
C 

12 CALL DIRNM (A» Mf D* T» MD) 

C ROOTS OF W INVERSE * A ARE IN T AND COLUMN EIGENVECTORS ARE IN D. 
C 

EM«M 

EKG=KG 

EN-N 

EKC=KC 

XL= 1 • 0 

TRACE«0.0 

DO 13 J * 1 *MD 

U( J)=T( J)/(1.0+T( JJ ) 

V( J>»SQRT<U<J> ) 

W( J) *1.0/ < 1.0+TIJ) ) 

XL=XL*W ( J ) 

13 TRACE=TRACE+T(J) 

D2= ( EN-lt-EKG*EM)*TRAC£ 

DO 14 J* 1 *MD 

14 2(J)=100«0*(T(J)/TRACE) 

IF ( M-2 > 16 • 16 * 17 

16 IF (XG-3) 1 6 » 18 ♦ 17 
18 YL=XL 

Fl*2»0 

F2*EN-3.0-EKC 
GO TO 19 

17 SL = SORT ( < (EM*EM)*nEKG“l«0>**2 ) -4«0 > / l ( EM*EM ) ♦ 

2 MEKG-1*0)**2)-5.Q) ) 

YL=XL** < 1 • 0 / SL ) 

PL« <EN-l«0-EKC)-{ (EM+EKG ) / 2 • 0 > 

QL= - ( ( EM * { EKG-1»0 ) )-2 *0 ) /4i0 
RL*tEM*<EKG-l«0) ) /2 *0 
Fl-2 • 0*RL 
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DISCRM CONTINUED 


F2=!PL*SL)+! 2 • 0*QL ) 

19 N 1 =F 1 
N2 a F2 

F*( <1.0-YLJ/YU*(F2/F1) 

YL* 1 • Q*»XL 

WRITE ! 6 # 201 ) XL # YL 

201 FORMAT ( 1 0W I LKS LAMBDA » '*F5*4#< GENERALIZED CORRELATION RATIO# 
CETA SQUARE = * #F5«4) 

WR I TE ( 6 » 20 ) F 

20 FORMAT 145H0F -RATIO FOR H2 ^OVERALL DISCRIMINATION#* F9.2) 

WRITE(6#21) Nl# N2 

21 FORMAT ( 8H0NDF 1 = 13# 12H AND NDF2 = 16) 

J=M D 

X! J+1)=1«0 

22 X( J)=X( J+l )*W! J) 

J*J-1 

I F t J ) 23»23#22 

23 DO 24 J* 1 * MD 

24 Y( J)= - PL * ALOG ( X ( J ) ) 

WRITE(6#1002) TRACE • D2 

1002 FORMAT !' QTRACE OF W-INVERSE*A= ♦•FlOtS//# 1 GENERALIZED MAHANANOBIS 
1 D-SQUARE * * #F15 • 5 ) 

WRITE ( 6 #25 ) 

25 FORMAT ( 4BH0CHI-SQUARE TESTS WITH SUCCESSIVE ROOTS REMOVED ) 
WRITE<6#261) 

261 FORMAT! 1HQ#20X*22H!ETA) (ETA SQUARE)) 

WR I TE ( 6 # 26 ) 

26 FORMAT! f OROOTS REMOVED CANONICAL R R SQUARED EIGENVALUE 
1CHI-SQUARE N#D*F. LAMBDA PERCENT TRACE ') 

DO 27 J = 1 #MD 
JT= J“1 

NDF= ( M-JT ) * ( KG-JT-1 *0 ) 

27 WRITE(6»28) JT*V< J) #U( J) #T( J) #Y( J) » NDF • X!j)#Z!j) 

28 FORMAT! 6X».I4#9X»2 (F6.3 »8X) #F9.3»5XfF10 t 0#4XfI5»2X»F9.2»F8#2) 



DlSCRM CONTINUED 


DO 29 J = 1 * MD 
DO 29 K,® 1 1 M 
A(J»KI=0.0 
DO 29 L = 1 *M 

29 A( J*K)-A( J*K)+D(L»J)*(C(L*K)/(EN-1.0) ) 

DO 30 J= 1 *MD 

DO 30 K a 1 *MD 
B(J*K)=0«0 
DO 30 L = 1 * M 

30 B( JtK)*BIJ»K)+A( JtL)*D<L*lO 
DO 31 J®1*M 

DO 31 (C = 1 *MD 

31 D( JtK)«D( J*K)*(1.0/SQRT(B<K*K)) ) 

WRITE(6»32) 

WR 1 TE ( 7 * 32 ) 

32 FORMAT ( 25HOROW COEFFICIENTS VECTORS ) 

DO 33 J B 1 f MD 

WR I TE 1 6 *49 } J* ( D ( K * J ) * K*1*MJ 

33 WRITE(7*49) Jf (D(KtJ) « K«1*M) 

49 FORMAT t 5H D F I3»2X*5E14*7 / ( 10X t 5E14.7 ) ) 

DO 34 J = 1 1 M 

34 2( J)*SORT(Ct JtJ) / (EN-1*0>) 

C TOTAL SAMPLE STANDARD DEVIATIONS ARE NOW IN 2. 

DO 35 J a 1 f M 
DO 35 K*1#M 

35 CtJ><)=C(Jf<) / ( £N*Z ( J ) *2 ( K. ) ) 

C TOTAL SAMPLE CORRELATION MATRIX IS NOW IN C. 

DO 36 J * 1 * M 
DO 36 K»1*MD 

36 B ( J*K ) »D ( J »K ) *Z t J ) 

DO 37 J» 1 1 M 

DO 37 K® 1 *MD 
A ( J«K ) ® 0*0 
DO 37 L = 1 * M 

37 At JtKI«AU*K )+C( J*L)*B tL*K ) 



U KJ 


D1SCRM CONTINUED 


WR I TE ( 6 * 38 ) 

38 FORMATS 42HOFACTOR PATTERN FOR DISCRIMINANT FUNCTIONS ) 

DO 39 J«1»M 

39 WR ITE ( 6 *40 ) J*!A(U*K>* K*1*MD) 

AO FORMAT ( 5H0TEST 1 4 # 1 Q ! 3X * F7 • 3 > / ( 9 X * 10 ( 3X *F 7 . 3 ) >) 

DO A 1 J*1*M 
T! J)*Q.O 
DO A1 K s 1 *MD 

A1 T( J)=T< J)+A( JtK)*A! J*K) 

WR I TE ! 6 *42 ) 

42 FORMAT! 19HOCOMMUNAL I T I TES FOR I5*21H DISCRIMINANT FACTORS ) 

WR I TE ( 6 1 43 ) (J* T ( J ) * J*1*M) 

43 FORMAT! 1H0*1Q!2X* 13* F7#3)) 

DO 44 J=1 »MD 

T(J)*0*0 
DO 44 K = 1 *M 

44 T!J) fi T{J)+A!K*J)*A(K»J) 

WR I TE ( 6 *45 ) 

45 FORMAT! 53H0PERCENTAGE OF TRACE OF R ACCOUNTED FOR BY EACH ROOT ) 
DO 46 J*1»MD 

46 T( J)*10Q*G#(T( J) / EM) 

WR I TE ! 6 » 43 ) !J* T C J ) * J*1*MD> 

KGT«KG+1 
DO 47 J = 1 »KGT 

47 RE AD (1*2301) ( A ( J »K ) *K«1 #M ) 

REWIND 1 

READS GROUP MEAN VECTORS AND GRAND MEAN VECTOR INTO COLUMNS OF At 
COLUMN KGT CONTAINS THE GRAND MEANS# 

DO 48 J*1 * KG 
DO 51 K “ 1 *MD 
T ( IO “0*0 
DO 51 L*1»M 

51 T!K)=T!K) + (A( J*L) - A ! KGT ♦ L ) ) * D(L*K) 

WR I TE ! 6 * 50 ) J* MD 
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DISCRM CONTINUED 


50 FORMAT (20H0CENTRQID FOR GROUP I 4*4H IN I4*32H DIMENSIONAL DISCRI 

2MINANT SPACE ) 

48 WRITE(6>43) (K* TOO* K=1*MD) 

NE = 1 

DO 1000 K 1 » 1 * KG 
DO 1001 J= 1 *M 

1001 RE AD (8*2301) ( A ( J * K ) *K= J * M ) 

DO 101 J=1 *M 
DO 101 K= 1 *M 

101 A ( K tU ) = A ( J *K ) 

DO 104 11=1 *MD 
DO 103 I 2 = 1 *M 
E( 11*12 )*0*0 
DO 102 I 3= 1 *M 

102 E(I1*I2)=E (11*12)^(0(13*11) *A (13*12)) 

103 CONTINUE 

104 CONTINUE 

DO 107 11=1 *MD 
DO 106 I 2 = 1 »MD 
A ( 11*12) =0*0 
DO 105 ! 3= 1 *M 

105 A(U*I2)»A(U*I2)+(E(U*I3)*D (13*12)) 

106 CONTINUE 

107 CONTINUE 
WRITE(6*2012)K1*IWS(K1) 

2012 FORMAT! /////IX* ‘DISPERSION IN REDUCED SPARE FOR GROUP ‘*12* 

1* WATERSHED **I5) 

DO 2013 1 1 = 1 *MD 

WRITE (6 >108) (A(U *12) , 12*1 *MD) 

108 FORMAT (10X*8(E14«7>2X) ) 

2013 CONTINUE 
1000 CONTINUE 

REWIND 8 

RETURN 

END 
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DIRNM 

SUBROUTINE DIRNM { A #M # B • X * XL ♦ LVECT ) 

C SUBROUTINE DIRNM# D I AGON AL I Z AT I ON OF A REAL NON-SYMMETR I C MATRIX 
C OF THE FORM B-INVERSE*A. CODED BY P* R* LOHNES# U# N • H# 

C A# M# B# X# AND XL ARE DUMMY NAMES AND MAY BE CHANGED IN THE 

C CALLING STATEMENT. 

C A AND B ARE M BY M INPUT MATRICES. UPON RETURN VECTOR XL CONTAINS 

C THE EIGENVALUES OF B"1#A# AND MATRIX X CONTAINS THE EIGENVECTORS 

C IN ITS COLUMNS. SUBROUTINE HOW PACKAGE IS REQUIRED. 

C LVECT SPECIFIES THE NUMBER OF EIGENVECTORS TO BE RETURNED. 

C 

DIMENSION A(20#20)#B<20»20)tX(20»20)fXL(20) #DUM1(2Q) #DUM2(20) 

1 # DUM3 12 0) ♦ DUMA (20) #E(2Q) 

C 

CALL HOW ( M# 20 #M # B #XL # X tDUMl #DUM2 #DUM3 #DUM4 } 

C 

DO 13 I«ltM 
DIAG*SORT(ABS(XLU J )) 

DO 13 J«1#M 
13 B ( J # I )=*X( Jtl >*DIAG 
16 DO 1 I * 1 #M 

1 XL C I ) «1 • O/SQRT ( ABS ( XL ( I ) ) ) 

DO 2 I » 1 #M 

DO 2 J«1#M 

2 B(I*J)«X( I • J ) *XL ( J ) 

DO 3 I * 1 #M 

DO 3 J=1 *M 
X ( I #J) = 0.0 
DO 3 K= 1 *M 

3 X( I #J)+B<K#I )*A (K * J J 

DO 4 I = 1 #M 

DO 4 J«1»M 
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DIRNM CONTINUED 



A { I f J ) =0*0 
DO 4 <= 1 #M 

4 A(ItJJ»A(I»J)+X(ItK)*B(K#J) 

C A NOW CONTAINS BP R I ME # A#B OF THE NOTES. 

C 

c 

CALL HOW(M»20*LVECT*A.XL*X*DUM1iDUM2 *DUM3*DUM4) 
C 

DO 6 I * 1 *M 
DO 6 J»ltM 
A ( I » J)»0«0 
DO 6 K= 1 »M 

6 A ( I # J ) »A ( I t J ) +B ( I *K)*X(KtJ) 

DO 9 I«1»M 

5UMV*0* 0 
DO 7 J» 1 »M 

7 SUMV s SUMV+ ( A U * I ) **2 ) 

DEN=SQRT (SUMV) 

DO 8 U«1»M 

8 XU. I J a A ( J » I ) /DEN 

9 CONTINUE 
C 

C COLUMNS OF XtZiJ) ARE NOW NORMALIZED# 

C 

RETURN 

END 
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HOW 


SUBROUTINE HOWIMVAR tMD IM *NVECT » R # E ♦ V • A # B * C * D ) 

EIGENVALUES AND EIGENVECTORS OF A SYMMETRIC MATRIX # ALGORITHM 
BY HOUSEHOLDER * ORTEGA * AND WILKINSON# ORIGINAL PROGRAM BY 
DAVID W# MATULA UNDER THE DIRECTION OF WILLIAM MERED I TH #UN I VERS I T Y 
OF CALIFORNIA# BERKELE Y # 1 96 2 . 

MODIFIED BY P. R. LOHNES# PROJECT TALENT# 1966# 

M IS THE ORDER OF THE INPUT MATR IX#R* 

MD IS THE DIMENSIONED SIZE OF R IN THE MAIN PROGRAM. 

NV IS THE NUMBER OF EIGENVECTORS TO BE COMPUTED# 

E IS THE VECTOR IN WHICH THE EIGENVALUES ARE RETURNED. 

V IS THE MATRIX IN WHICH THE EIGENVECTORS ARE RETURNED. 

THE EIGENVECTORS ARE STORED AS COLUMNS IN V. 

A #9 >C #AND D ARE WORKSPACE VECTORS. 

DIMENSION R( 1 ) *E( 1 ) *V( 1 ) #A(1 ) #B( 1 1 #C ( 1 ) #D( 1 ) 

C 

M=MVAR 

MD-MDIM 

NV=NVECT 

IF(M-l) 100#97»96 
96 

C TRI-DI AGONAL IZE THE MATRIX. 

M2*M1*MD+M 

M3*M2-MD 

MA=MD+1 

L*0 

DO 1 1*1 *M2 #M4 
L“L+1 

1 A ( L ) =R ( I ) 

B ( 1 ) =0. 0 
I F ( M“2 ) 1 3 * 2 # 3 
3 KK = 0 

DO 15 K a 2 #M1 
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HOW CONTINUED 



K.L = KK+K 
KU=KK+M 
KJ«K+1 
SUM= 0.0 
DO 4 J»KL *KU 

4 SUM=SUM+R( J)**2 
S*SQRT( SUM) 

Z=R(KL5 

B { K ) «S I GN ( S *-Z ) 

S=1.0/S 

CUO=SQRT(ABS (Z) * S + 1.0) 
X»SIGN(S / C(K) * Z) 

RUL)*C(K) 

DO 5 I«K J *M 
J J= I +KK 
C( I >»X*R( JJ) 

5 RIJJK(I) 

DO 8 J = K*M 
JJ*J+1 

D ( J) “0. 0 
L-KK+J 
DO 6 1*K*J 
L=L+MD 

6 Dl J)=D( J)+R(L)#C< I ) 

IF(JJ-M) 7*7*9 

7 DO 8 1 « J J *M 
L«L + 1 

8 D( J)=D( J)+R(L)*C( I) 

9 X = 0.0 

DO 10 J*K*M 

10 X«X+C(J)*D{J) 

X« • 50#X 

DO 11 I =K * M 

11 D< I ) *X#C < I )-D(l ) 

LL*KK 
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HOW CONTINUED 


KK=KK+MD 
DO 15 I*K#M 
LL=IL+MD 
DO 15 J=I*M 
L”LL+J 



15 

R(L)=R(L)+D( I )*Cl J)+Dl J)*C< I ) 



L“1 

DO 12 I » 1 * M 
X«A ( I ) 

o o 


A ( I ) »R C L > 



R ( L ) a X 

►a 3 

12 

L* L+M4 

gs 

2 

B ( M ) »R { M3 > 


C 

COMPUTE EIGENVALUES* 


13 

BD=ABS (All ) ) 

g£ 


DO 14 I «2 >M 

fa § 


IF (BD-IABSIAU ) )+B(I)^2 ) ) 130 #14# 14 

130 

BD= ABS ( A ( I > ) +B ( I )**2 

« m 

14 

CONTINUE 


BD=BD+i • 0 


DO 1.6 I*1#M 
A(I)«A(I)/BD 

b< n«B( n /bd 

D(I)«liO 

16 El I ) *-l *0 
DO 37 K«1 »M 

17 IF (ABSlDlK) J-ABSlElK) ) )101*101»105 
101 AMAX1= ABSIEIK) } 

GO TO 110 

105 AMAX1* ABS(DIK) ) 

110 IF IAMAX1-1.E-9) 115tl20«120 
115 AMAX1 = 1 • E-9 
120 CONTINUE 

IF I(DtK)-EIK) )/AMAXl-l*E"6) 37t37*18 

18 Xs(DlK)+EU) )*.5Q 


A58 



HOW CONTINUED 


I S2 = l 
S2= I S2 
C ( 1 ) *A( 1 )“X 
IF(CU) ) 19 *20»2Q 

19 IS1=-1 
S1*IS1 
N»0 

GO TO 21 

20 IS1-1 
S1-IS1 
N«1 

21 DO 31 I»2tM 
IF(B( I) ) 22 >26*22 

22 IF(Btl-l) )23»27*23 

23 IF ( ABS ( C ( 1-1 ) )+AB5lC( 1-2 ) )-l .OE-15 ) 24*25 * 25 

24 Cn-X)*C(I-l)*l*0E15 
C ( 1-2 ) *C ( I~2)*1*0E15 

25 C(I) ■ (A( I I-X)*CII“1)“B( I )**2 * C(I-2) 

GO TO 28 

26 C(I) = (A(I)-X) * SIGN CleO* SI) 

GO TO 28 

27 C(I) * { A ( I ) -X ) * CII-1) - SIGN (B(I)**2# S2 ) 

28 S2 * SI 
I S 2 = S 2 

IF (Cdl) 29* 30# 29 

29 SI « SIGN (S1*C( I) ) 

IS1=S1 

IF (IS2 + I5U 30# 31# 30 

30 N = N + I 

31 CONTINUE 
N = M - N 

IF IN - K> 34# 32# 32 

32 DO 33 J = K. N 

33 D ( J ) * X 

34 N = N + 1 
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HOW CONTINUED 


IF ( M - N) 1 7 1 35* 35 

35 DO 36 J * N* M 

IF (X - £(J}) 17* 17* 36 

36 EU) = X 
GO TO 17 

37 CONTINUE 

DO 38 I = 1# M 
All) « Ail) * BO 
811) * BU) * BD 

38 CII) ■ (DID + EU)) * BD * .50 
Ml = M 

K = 1 

39 I * 1 

40 DO 43 J = It Ml 

IF (I - J) 41* 43* 41 

41 IF (CII) "C ( J ) ) 43* 43* 42 

42 I - J 

GO TO 40 

43 CONTINUE 

eu> a cm 

K = K + 1 
Ml = Ml - 1 

IF (I - Ml - 1) 44* 46* 46 

44 DO 45 M2 = 1* Ml 

45 C ( M2 ) = CIM2+1) 

46 IF (Ml - 1) 47* 47* 39 

47 EU) * Cil) 

IF (ISIGN (1* NV)) 79* 76* 76 

76 DO 77 I * 1 * M 

77 C(l) « Eli) 

J = M 

DO 78 I = 1* M 
Ell) * C( J) 

78 J s J " 1 

79 CONTINUE 
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HOW CONTINUED 


DECIDE WHETHER TO COMPUTE EIGENVECTORS* AND IF SO* HOW MANY 

IF (NV) 48* 99* 48 

48 KX = IABS(NV) 

J = 1 

DO 98 INV * l*KX 
X * A ( 1 ) - E (INV) 

Y a B ( 2 ) 

Ml a M - 1 

DO 54 I = 1* Ml 
IJ = J + I « 1 

IF ( ABS (X) - ABS (B(I+1))> 49* 51* 53 

49 C( I) = B( 1 + 1 ) 

D ( I ) = A ( I + 1 ) - EUNV) 

V( IJ) = Bt 1 + 2) 

Z = -X / CII) 

X=Z*D(I)+Y 

IF (Ml - I) 50* 54* 50 

50 Y = Z * V( IJ) 

GO TO 54 

51 IF (X) 53* 52* 53 

52 X - 1.0E-10 

53 C( I ) = X 
D( I) * Y 
V(IJ) * 0.0 

X = A ( I + 1 ) - IBII+I) / X * Y + EUNV)) 

Y = B ( 1 + 2 ) 

54 CONTINUE 

MJ n M + J - 1 

IF (X) 56* 60* 56 

56 V ( MJ ) = 1.0 / X 

57 I = Ml 

IJ a J + I - 1 

VUJ) a (1.0 - D ( I ) * V ( M J ) ) / C(I) 
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HOW CONTINUED 


X = V<MJ)**2 + VI I J )**2 

58 r = i - i 

IJ ■ J + I - 1 

IF (I) 59 * 6 1 f 59 

59 VUJ) « (1*0 - (D(I) * VUJ+1J + VUJ) * V ( I J+ 2 ) ) ) / C(I) 
X = X + V( I J ) **2 

GO TO 58 

60 V(MJ) = ItOElO 
GO TO 57 

61 X 8 SORT ( X ) 

DO 62 I * It M 
IJ ■ J + I - 1 

62 VUJ) ■ VIIJ) / X 
J1 = Ml * MD - MD 
K = M 

GO TO 66 

63 K 8 K - 1 

Jl * J1 - MD 
Y = 0*0 

DO 64 I « Kt M 
IJ * J + I - 1 
L = Jl + I 

64 Y = Y + V( IJ) * R(L) 

DO 65 I * K * M 

IJ = J + I - 1 
L ! Jl + I 

65 VUJ) * VUJ) - Y * R(L) 

66 IF (Jl) 63* 67 » 63 

67 NPIUS • 0 
NMIN * 0 

DO 70 I ■ 1# M 

IJ * J + I - 1 

IF (V(IJ)) 68* 69* 69 

68 NMIN * NMIN ♦ 1 
GO TO 70 
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HOW CONTINUED 



69 NPLUS * NPLUS + 1 

70 CONTINUE 

IF (NPLUS - NM IN ) 71* 73* 73 

71 DO 72 I = 1* M 
1J s J + I * 1 

72 VII J) * -V ( I J) 

73 CONTINUE 

98 J = J + MO 

C RESTORE THE INPUT MATRIX. 

99 MD1 = MD + 1 
JJ = MD1 

Ml » M * MD 

DO 75 I * 2* Ml* MD1 

K = I 

DO 74 J = JJ* Ml* MD 
R(K) = RIJ) 

74 K = K + 1 

75 JJ - JJ + MD 1 
GO TO 100 

97 Ed) = R(l) 

V(l) = 1.0 
100 RETURN 
END 
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MANOVA 


C MULTIVARIATE ANALYSIS OF VARIANCE# A CQGLEY-LQHNE S ROUTINE# 

C 

C THIS PROGRAM COMPUTES MANOVA TESTS OF HI { EQUALITY OF DISPERSION 

C AND H2 (EQUALITY OF CENTROIDS)# UNIVARIATE F-RATIQS FOR MEANS# 

C SELECTED SAMPLE STATISTICS* AND THE W (POOLED WITHIN-GRQUP SSCP) 

C AND T(TOTAL SAMPLE SSCP) MATRICES REQUIRED FOR THE DISCRIMINANT 
C ANALYSIS PROGRAM# THESE MATRICES ARE PUNCHED IN UPPER-TRIANGULAR 

C FORM# THE PROGRAM WILL PROCESS UP TO 50 VARIABLES AND ANY NUMBER 

C OF GROUPS* 

C 

C INPUT 

C 

C 1 ) FIRST TEN CARDS OF THE DATA DECK DESCRIBE THE PROBLEM IN A TEXT 

C WHICH WILL BE REPRODUCED ON THE OUTPUT. DO NOT USE COLUMN 1. 

C 2) CONTROL CARD ( CARD 11) 

C COLS 1-2 M*NUMBER OF VARIABLES 

C COLS 3-5 KG=NUMBER OF GROUPS 

C 3 ) FORMAT CARD ( CARD 12) 

C 4 ) EACH GROUP OF SCORE CARDS IS PRECEDED BY A CARD GIVING 

C NG=NUMBER OF SUBJECTS IN THE GRQUP(COLS 1-5)# 

C THUS# SUBJECTS MUST BE SORTED INTO GROUPS AND THE GROUPS COUNTED 

C BEFORE MANOVA CAN BE RUN# 

C 

C PUNCHED OUTPUT IS ALL TO FORMAT ( 10X * 5E14# 7 /(10X# 5E14.7M.AND IS 

C 

C 1 ) GROUP MEANS# FOLLOWED BY GRAND MEANS# 

C 2 ) T MATRIX(T0TAL SAMPLE DEVIATION SSCP MATRIX) 

C 3)W MATRIX (POOLED W I TH IN-GROUPS DEVIATION SSCP MATRIX) 

C 4)D INVERSE( INVERSE OF POOLED-SAMPLES DISPERSION t ST i M aTl / 

C 

C SUBROUTINE MAT I NV IS REQUIRED. 

C 

C 

DIMENSION T IT ( 10*20) #A( 20 #20 ) #B( 20.20 ) *C( 20.20 J # 
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MANOVA CONTINUED 


*T ( 20 ) *U ( 20 ) *W ( 20 ) t X l 20 ) *0120*20) #NV ( 20 ) * 

* Z WS C 20) * ISCN( 20) *SDA(4*4) * SM ( 4 ) * Y ( 3 ) 

COMMON M*KG*N»KC*TITtNTrTLE*IWS» ISCNtNV 
READ(5tll21) I GROUP 
X 1 2 1 FORMAT ( 12) 

C 

DO 1122 1 = 1*1 CROUP 
RE AD (5*1101) I SN * I SHD ♦ I T I L *NG 
1101 FORMAT (6X*I5*T19*I5*T31»I5*T42*I8) 

Y ( 1 ) = I SN 
Y ( 2 ) * I SHD 
Y ( 3 ) =NG 
WRITE(2*102)Y 
DO 1123 J= 1 * A 

READ (5*1102) ( SDA( J*K) *K=1 *4 ) 

WRITE (2 *102) <$DA( J#K) *K=lt4) 

1123 CONTINUE 

RE AD (5*1102) ( SM ( K ) *K*1 *4 ) 
WRITE(2*102MSM(K)*K-1*4> 

1122 CONTINUE 
C 

REWIND 2 
C 

1 WR ITE (6*2) 

2 FORMAT ( 3 3H1 MANOVA • A COOLE Y-LOHN ES PROGRAM 
READ ( 5* 1001 * END = 53 ) NT I TIE 

DO 3 J= 1 *NT I TL£ 

READ (5*4) (TIT(J*<)*K=1*20) 

3 WRITE(6*4) (TlT( J*K) »K=1*20) 

4 FORMAT ( 20A4 ) 

READ (5*5) M * KG 

5 FORMAT ( 12*13) 

READ (5*1001) ( NV ( I ) *I=1*M) 

1001 FORMAT (2513) 

WR I TE ( 6 » 1003 ) ( NV ( I ) *I = 1*M) 



MANOVA CONTINUED 


c o 


sg J“ 

§ 
E3 


1003 FORMAT! ‘OVARIAaLES FOR THIS RUN ARE * • / * 2 5 I 3/ * 2 8X * 2 5 1 3 ) 

102 FORMAT ( 162 ! 10 A4 ) ) 

1102 FORMAT! 19X*4E14*7) 

EM = M 
E1CG*KG 
EK = KG 

WRITE!6i6) M#KG 

6 FORMAT! 13H0ANALYSIS FOR I3*14H VARIABLES ANDI4*7H GROUPS) 
WR I TE ( 6 *9 ) 

9 FORMAT ( 1H0 • 25 ( 5H - ) ) 

DO 7 J-1*M 
T ! J ) =0* 0 

■ DO 7 K= 1 *M 
B!J*K)=0*Q 

7 C! J»K)=0.0 
H1LOG$»0 • 0 
GA1S = 0* 0 
FA1S B 0#0 
N*0 

c 

DO 19 I G=1 »KG 
READ ! 2# 102 ) Y 
I SCN ( I G ) S Y ( 1 ) 

IWS! 1G) =Y( 2 ) 

NG=Y( 3) 

DO 1124 1-1*4 

1124 READ! 2*102 ) ! SDA! I .<) »K=1 *4) 

READ (2*102) { SM ( K ) *K=1*4) 

ENG*NG 
N»N+NG 
WRITE (6*9) 

WRITE<6#10) IGtNG 

10 FORMAT! /1X» ’GROUP 1 *13* 1 NG*‘*I7) 

WRITE(6« 1104) ISCN! IG) * IWS( IG) 

1104 FORMAT! /IX* » ERTS SCENE SIS* 1 WATERSHED »* I 5 ) 
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ORIdNAU PAGE IS 

of poor guAura 


MANOVA CONTINUED 


DO 11 J»1#M 
Ut J)»SM(NV< J) ) 

DO 11 K«1*M 

11 A( JiK)=SDA(NV( J> *NV(<) ) 

DO 12 J*1»M 
TIJ)»T(J)+U(J) 

DO 12 K*1#M 

12 C( J*K)=C( JtK)+A( JtK) 

DO 13 J ® 1 • M 

DO 13 K = 1 * M 

A ! J*K ) sA ( J*K J-Ut J ) * U(K) / ENG 
B( J*tO=B( J#K)+A( J.K) 

13 A( J#K)=A( J#K ) / ! ENG-1 • 0 5 
DO 14 J = 1 * M 

U ( J ) S U ( J ) /ENG 

14 W( J)*SQRT(A! JtJ5 ) 

WRITE!6 >15 ) IG 

13 FORMAT! 17H0MEANS FOR GROUP 14) 

WRITE(6*16) ( U(JJt J = 1 > M) 

16 FORMAT ( 1H0*1Q(3X»F7.2) ) 

WRITEC1 » 1 02 ) (U! J) »J = 1»M) 

30 FORMAT { 4H R0WI3|3X* 5E14.7 / ( 1 OX * 5E 14 . 7 ) ) 
WRITE!6tl7) 

17 FORMAT! 21H0STANDARD DEVIATIONS ) 

WR I TE ( 6 * 1 6 ) ( W ( J ) .J*ltM> 

DO 1004 J= 1 # M 

1004 WRITE!8»102) ! A ! J * K, ) *K=J*M) 

CALL MATINV! A»M*DET ) 

WR I TE ( 6 ► 1 8 ) DET 

16 FORMAT { 26HODI5PERSION DETERMINANT « E14.4) 

H1LOGS-H1LOGS+! ! ENG-1.0) *ALOG { DET ) ) 
FA1S=FA1S+! 1.0/! ENG-1 .0 ) ) 

GAlS a GAlS+(1.0/! ! ENG-1.0 J*#2 ) ) 

19 WR I TE ( 6 1 9 ) 
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KANOVA CONTINUED 



REWIND 2 
REWIND 8 
C 

c 

EN S N 

DO 20 J = 1 1 M 
DO 20 K*1 *M 

A( J*K)=C( J#iO“T< J) * T(K)/EN 
D { J#K ) «A ( J * K ) 

20 C( J*K)*B( J*K) / (EN-EKG) 

DO 21 J 3 1 * M 

T( J)»T< J)/EN 

21 U( J)«SQRT (C(JiJ) ) 

WR I TE { 6 ♦ 22 ) 

22 FORMAT ( 23H0MEANS FOR TOTAL SAMPLE) 

WR I TE { 6 * 1 6 ) ( T( J) • J=1»M) 

KGT=KG + 1 

WR I TE ( 1 * 102 ) { T { J } tJ*ltM) 

REWIND 1 
WR ITE (6 * 23 ) 

23 FORMAT { 35HOPOOLED-5AMPLES STANDARD DEVIATIONS ) 

WRITEJ6 f 16 HU(J) * J*ltM5 

WR I TE { 6 * 9 ) 

WRITE<6*38) 

38 FORMAT ( 9H0T MATRIX) 

DO 34 J= 1 »M 

WRITE (6 *30) J* (A( J»K) *K*J»M) 

34 WRITE(3*102) (A( J*K) tKsJ.M) 

REWIND 3 

WR I TE ( 6 * 9 ) 

DO 35 J * 1 *M 
DO 35 K = 1 *M 

35 A* J*K)*A( J#K) 

C A IS NOW THE A ( AMONG-GROUP S SSCP) MATRIX* B IS NOW THE w (wIThIn 

C GROUPS SSCP) MATRIX. C IS NOW THE POOLED-GROUPS DISPERSION tST. 
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MANOVA CONTINUED 


C 

WRITE(6*28 ) 

28 FORMAT ( 9H0A MATRIX) 

DO 29 J= 1 *M 

29 WR I TE ( 6 * 30 ) J* (A(JiK)t K*J*M> 

WR I TE { 6 * 9 ) 

C 

WRITE(6»31) 

31 FORMAT { 9H0W MATRIX) 

00 32 J = 1 * M 

WRITE (6 *30) J» (B( J*K) *K=J*M> 

32 WRITEU.102) <B( J*K) .K=J.M> 

REWIND A 

WR I TE ( 6 * 9 ) 

C 

CALL MATINV(C«M»DET) 

C DISPERSION MATRIX STORED ON DISK FILE. 

WR I TE ( 6 * 1006 ) 

1006 FORMAT ( * OC MATRIX*) 

DO 33 J=1»M 

33 WRITE(6*30)J*(C(J»K)«K=J*M) 

C 

WR!TE(6*18) DET 

HlLOG=( EN-EK) * ALOG (DET) 

XMM*H1L0G-H1L0GS 

FI®. 5 * (EK-1.0) * EM * (EM+1.0) 

A1A=( FAlS-( 1.0/(EN-£K) ) )*l(2i0 * ( EM*EM ) ) + ( 3 • 0 * EM)-1.0) 

A 1 «A1A/ ( 6 • 0 * (EK - 1.0) * (EM+l.QM 

A2® { GA1S -(1.0 / ( EN-EK ) **2 ) ) * ((EM-1.0) * (EM+2.0)) 

C/ ( 6 • 0 *( EK-1.0) ) 

DIF = A2-A1 * A 1 
IF(DIF) 24*24*25 
24 P2*(Fl+2.0)/(Al * A1-A2) 

B1«F2 / (1.0 -Al+(2.0/F2) ) 

F»(F2 * XMM ) /(FI * (81 - XMM ) ) 
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MANOVA CONTINUED 


C 


c 

c 

c 


GO TO 45 

25 F2*(F1+2#Q)/DIF 

B1*F1 / ( 1 • Q“A1 - (Fl / F2 ) ) 


F«XMM/B1 
45 NDF1»FX 
NDF2 ®F2 

WR I T£ ( 6 * 26 ) XMM*F 

26 FORMATS 47H0F0R TEST OF Hi 

C10H AND F * F10.3) 

WR ITE ( 6 *27 ) NDF1 1 NDF2 

27 FORMAT ( 15HOFOR F* NDFi « 
WR I TE ( 6 »9 ) 


( equal i t v of dispersions)# m 

3* 12H AND NDF2 * 19) 


Fl 0 • 3 « 


N1*EKG-1*Q 
N2*EN-EKG 
WR I TE ( 6 * 9 ) 

WRITE(6*4Q)N1»N2 

40 WRi?El^9y° UNIVARIATE F_RATI0S ’ WITH « DF1 = I3.12H AND NDF2 = 16) 
WRITE<6»41) 

41 FORMAT) 71HOVARIABLE AMONG MEAN SO WITHIN MEAN SQ F-RATIO 
C ETA SQUARE) 

DO 42 J = 1 * M 

ET ASQ s A ( J * J ) / ( A ( J * J ) + 9 ( J » J ) j 
AMS*A { J # J ) / ( EKG-1 »0 ) 

WMS = B f J # J ) / (EN-EKG) 

F« AMS/WMS 

42 WRITE (6 *43 ) J * AMS * WMS *F *ETA5G 

43 R °^*|' 3 ^* I3,5X » F 9*2.UX,F9.2.10X.F7.2.8X.F5.4) 

CALL MATINV( B*M*DETW) 

CA^L MAT?N^DyM^E^T? F POOLED - SAMPLES DEV I AT I ON SSCP MATRIX, W. 
DETT IS DETERMINANT OF TOTAL SAMPLE DEVIATION SSCP MATRIX* T» 
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MANOVA CONTINUED 


XL-DETW/DETT 
YL S 1 • 0-XL 
WRITE(6*46)XL*YL 

46 FORMATS * OWI LKS LAMBDA = ' #F7 .4* 1 GENERALIZED CORRELATION RATIO* 

1 ETA SQUARE = 1 *F5.4) 

I F ( M“2 ) 47 * 47 *49 

47 IFCKG-3) 48*48*49 

48 YL=XL 
Fl*2 • 0 

F2 = EN-3 #0 
GO TO 50 

49 SL=SORT ( ( ( EM * EM) * ( ( EKG - 1.0)**2) - 4*0) / UEM * EM) + 

2 I (EKG - l.Q)**2) - 5*0) ) 

YL*XL ** (1*0 / SL J 

PL S ( EN-1 *0 ) - ( ( EM+EKG ) / 2.0) 

QL*-((EM * ( EKG -l.Q))-2.Q) / 4.0 
RL“ ( EM * (EKG-1.0) ) / 2.0 
Fl=2.0 * RL 

F2*(PL * SU + 12.0 * QL) 

50 N1*F1 
N2 = F2 

Fa(d.O-YL) / YL) * IF 2 / FI) 

WR ITE ( 6 *51 ) F 

51 FORMAT(45HOF-RATIO FOR H2 » OVERALL DISCRIMINATION* = F9.2) 

WR I TE ( 6 * 5 2 ) N 1 * N 2 

52 FORMAT { * 0NDF1 = dlS* 1 AND NDF2 * ‘ *19) 

WR I TE ( 6 * 9 ) 

KC S Q 

CALL DISCMX 
GO TO 1 

53 STOP 
END 
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MAT I NV 


*I-W 



SUBROUTINE MAT I NV ( A *M * OET ) 

DIMENSION AC 20 *20) *IPVT(100) *PVT (100 ) * IND C 100 * 2 ) 
DET«1« 

DO 1 J= 1 *M 

1 1 PVT ( J ) ®0 
DO 10 I « 1 #M 
AMAX=0 • 0 

DO 5 J* 1 *M 

IF( IPVTC J)-l ) 2*5*2 

2 DO 5 K=1*M 

I F ( IP VT (0-1)3*5*20 

3 IF ( ABS( AM AX ) "ABS ( A ( J iK ) ) )4*5*5 
A I ROW® J 

I COL=K 
AMAX*A ( U *K ) 

5 CONTINUE 

I PVT ( ICOL) = IPVT( !CQL)+1 
IF I IROW-ICOU6*8t6 

6 DET*-DET 
DO 7 L® 1 *M 
SWAP® A ( I ROW * L ) 

A ( IR0W#L ) «A ( I COL * L ) 

7 A ( I COL * L ) ® SWAP 

8 1 ND ( I # 1 ) * I ROW 
I ND ( I #2 ) * I COL 

PVT( I)«A( ICOL*ICOL) 

DET*DET*PVT( I) 

A ( I COL* I COL ) = 1 • 

DO 9 L“ 1 *M 

9 A( ICOL*L}=A( ICOL*L)/PVT( I ) 

DO 10 L1=1*M 

IF (Ll—ICOL) 11*10*11 
11 SWAP*A ( L 1 * I COL ) 

A ( LI * ICOL ) =0*0 
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MAT I NV CONTINUED 


DO 12 L = 1 «M 

12 A(Ll»U=A(LhL)-A( ICOL#L)*SWAP 
10 CONTINUE 

DO 20 I * 1 * M 
L=M+1-I 

IF(IND(L»1)-IND(L»2))13*20*13 

13 I ROW* IND C L * 1 ) 

I COL* IND ( L * 2 ) 

DO 20 K = 1 »M 
SWAP* A l K * I ROW ) 

AU*IROW)*AUf ICOL) 

A ( Kt I COL ) =SWAP 

20 CONTINUE 
RETURN 
END 
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SPLIT 


ENT 

SPLIT DC 
STX 
STX 
STX 
LDX 
LD 
STO 
LD 
STO 
LD 
SLA 
STO 
STO 
MDX 
STX 

ADD1 LDX 
ADD2 LDX, 
COUNT LDX 
LOOP SLT 
LD 
RTE 
SLT 
STO 
MDX 
SLA 
SLT 
STO 
MDX 
MDX 
MDX 
MDX 

XR1 LDX 
XR2 LDX 


SPLIT 

*-# 

1 XR1 + 1 

2 XR2 + 1 

3 XR3+1 
II SPLIT 

X 0 

ADD1 +1 
1 1 

ADD2+1 

112 

1 

II 2 

COUNT+1 
1 3 

1 BACK+1 
LI *-* 

L2 *-* 

L3 *“* 

16 
1 0 
16 
8 

2 0 
2 -1 

16 

8 

2 0 
1 -1 
2 -1 
3 -2 
LOOP 
LI *-* 

L2 
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SPLIT CONTINUED 


XRO 

BACK 


LDX 

BSC 

END 


L 3 *-# 
L *-* 
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HALF 


ENT 

HALF DC 
STX 
. STX 
STX 
LDX 
LD 
STO 
LD 
STO 
LD 
STO 
SLA 
STO 
MDX 
STX 

ADD1 LDX 
ADD2 LDX 
COUNT LDX 
LOOP LD 
SRT 
STO 
XCH 
SRT 
STO 
MDX 
MDX 
MDX 
MDX 

XR1 LDX 
XR2 LDX 
XR3 LDX 
BACK BSC 
END 


HALF 


1 

XR 1 + 1 

2 

XR2 + 1 

3 

XR3 + 1 

a 

HALF 

i 

0 


ADD1+1 

l 

1 


ADD2+1 

a 

2 


COUNT+1 


1 

a 

2 

l 

3 

i 

BACK+1 

Ll 


L2 


L3 

*-* 

1 

0 


€ 

2 

0 


8 

2 

-1 

1 

-1 

2 

-2 

3 

-1 


LOOP 

Ll 

#-* 

L2 

#-# 

L3 

*-* 

L 
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